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INTRODUCTION 


This  report  summarizes  the  results  of  a  study  to  assess  the  effects 
of  turbulent  fluctuations  on  the  infrared  thermal  emission  from  a  rocket 
exhaust  plume.  The  propellant  composition,  length  scales,  and  flight 
conditions  are  typical  of  those  for  a  small  tactical  missile.  We 
consider  emission  from  the  COg  molecule  in  the  4.3  pm  wavelength  region. 
This  chemical  species  is  formed  in  the  afterburning  region  where  the 
high  velocity,  high  temperature,  fuel  rich  rocket  exhaust  mixes  with  the 
ambient  air.  The  region  where  the  plume  and  ambient  air  mix  is  highly 
turbulent,  with  random  fluctuations  in  temperature  and  radiating  species 
concentration.  The  goal  of  this  study  was  to  determine  how  these 
fluctuating  properties  affect  the  gaseous  thermal  emission  from  this 
class  of  flows. 


Our  approach  was  to  utilize  the  available  turbulent  flowfield 
computational  capability  that  has  evolved  from  the  higher-order 
turbulence  modeling  efforts  at  Aeronautical  Research  Associates  of 
Princeton. 1-3  The  computational  codes  developed  during  these  studies 
predict  the  second-order  correlations  of  the  fluctuating  temperature  and 
chemical  species  in  addition  to  the  mean  (time-averaged)  properties  of 
the  flow.  These  predictions  allow  the  time-averaged  radiance  to  be 
evaluated  with  contributions  from  these  fluctuations  as  well  as  the  mean 
properties  of  the  exhaust  plume  flow. 

The  results  of  the  study  are  presented  in  the  following  sequence.  A 
brief  background  and  motivation  for  the  work  is  provided  in  the 
remainder  of  the  Introduction,  The  second  section  presents  a  discussion 
of  the  flowfield  calculations,  which  includes  a  systematic  comparison  of 
the  higher-order  closure  results  with  those  obtained  from  the  BOAT 


‘Donaldson,  C.  duP.,  "Atmospheric  Turbulence  and  Dispersal  of 
Atmospheric  Pollutants,"  Environmental  Protection  Agency,  March  1973. 
(EPA-R4-73-0 16a,  publication  UNCLASSIFIED.) 

2Fishburne,E.S. ,  R.A.Beddini,  and  A.K.  Varraa,  "The  Computation  of 
Afterburning  Rocket  Exhaust  Plumes,"  Aeronautical  Research  Associates  of 
Princeton,  Inc.,  July  1976.  (A.R.A.P.  Report  No.  283,  publication 

UNCLASSIFIED.) 

3Varma,  A.K.,  E.S.  Fishburne,  and  R.A.  Beddini,  "A  Second-Order  Closure 
Analysis  of  Turbulent  Diffi'ion  Flames,"  NASA  Langley  Research  Center, 
June  1977.  (NASA  CR-1 4522b,  publication  UNCLASSIFIED.) 
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module  of  the  JANNAF  Standard  Plume  Flowfield  code  (SPF). 4  This 
comparison  is  included  because  the  SPF  will  be  the  industry-wide  tool 
for  predicting  exhaust  plume  flows.  The  SPF  predictions  provide  a 
reference  against  which  the  more  detailed,  second-order  closure  results 
can  be  compared.  We  also  include  in  this  section  an  assessment  of  the 
importance  of  the  turbulent  fluctuations  on  the  chemical  reactions.  A 
derivation  of  the  equation  for  the  fluctuation  augmented  radiance  is 
given  in  the  third  section.  The  results  of  the  radiance  predictions  are 
given  and  discussed  in  the  fourth  section.  In  all  cases,  the  radiance 
with  fluctuations  is  compared  with  that  ignoring  the  fluctuations.  A 
summary  and  the  conclusions  are  given  in  the  final  section.5 


BACKGROUND 

Infrared  signature  codes  currently  available  for  routine  use, 
now5  or  in  the  near  future,6  account  only  for  emission  from  the  mean 
properties  of  the  flow.  However,  the  highly  non-linear  dependence  of 
the  emission  and  absorption  on  temperature  suggests  that  the  true, 
time-averaged*  radiance  should  be  different  from  that  evaluated  at  the 
mean  properties: 


N  (T,c)  i  N  (T,c)  (1) 

w  u 


A  simple  quantitative  example  is  easy  to  construct.  Assume  that  c.he 
temperature  dependence  of  the  Planck  function  BW(T)  can  be  approximated 
by  a  simple  power  law  of  the  form  aT10  (m  =  4  exactly  for  the  total 
radiance,  and  2  <  m  <  10  for  the  spectral  radiance,  depending  on  the 
wavelength  region  and  temperature).  The  time-averaged  radiance, 

4Da3h, S.M.,  et  al.,  "Operational  Instructions  for  a  Preliminary  Version 
of  the  JANNAF  Standard  Plume  Flowfield  Model  (SPF/1),"  June  1980. 
(A.R.A.P.  Report  No.  415,  publication  UNCLASSIFIED.) 

5Pergament,  H.S.,  et  al.,  "The  Naval  Weapons  Center  Target  Signature 
Code,"  February  1979.  (A.R.A.P.  Report  No.  380,  publication 

UNCLASSIFIED.) 

6Ludwig,  C.  et  al.,  "Development  and  Validation  of  a  Standardized 
Infrared  Radiation  Model  (SIRRM) ,"  JANNAF  Standard  Plume  Model  Workshop, 
U.S.  Army  Missile  Command,  Huntsville,  AL,  2-3  April  1981.  (paper, 
UNCLASSIFIED.) 

•Averaged  over  times  long  compared  to  the  time  scale  of  the  turbulent 
fluctuations. 
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acsounting  for  the  fluctuations  T*  about  the  mean  temperature  T  is, 
approximately. 


N  (T)/N  (T)  =  1  +  1/2  m(m  -  1)  Hi  +  +  ...  .  (2) 

w  w  J2  f4 


For  afterburning  rocket  plumes,  and  many  other  turbulent  flows, 
TT2/T2  ~  0.1,  suggesting  that  the  temperature  fluctuations  alone  can 
substantially  enhance  the  radiance  over  that  evaluated  at  the  mean 
temperature.  For  example,  if  m  =  5,  the  spectral  mean  radiance  i3 
doubled  by  the  temperature  fluctuations.  Of  course,  the  actual 
situation  in  a  partially  transparent  volume  of  molecular  emitters  is 
more  complicated.  The  species  and  temperature  sensitivity  of  the 
molecular  absorption  can  be  either  positive  or  negative,  and  there  are 
spatial  correlations  between  fluctuations  along  the  line  of  sight  that 
both  augment  and  diminish  the  radiance.  A  precise  evaluation  of  the  net 
change  in  radiance  requires  that  these  additional  factors  be  included. 
That  is  the  goal  of  this  study. 


FLOWFIELD  CALCULATIONS 


REVIEW  OF  AVAILABLE  CODES 

Two  state-of-the-art  computer  programs  for  the  calculation  of 
turbulent  reacting  flows  are  available  to  us  for  predicting  the 
flowfield  properties  of  missile  plumes  of  interest.  The  studies  are 
restricted  to  flowfields  involving  constant  pressure  mixing  and 
reactions.  One  is  the  BOAT7  code  -  the  constant  pressure  turbulent 
mixing  module  of  the  JANNAF  Standard  Plume  Model.  The  other  code  is  the 
A.R.A.P.  reacting  shear  layer  (RSL)  code.2  The  program  has  also  been 
called  the  multi-equation  turbulent  £eacting  code  (METREC)  in  some 
previous  publications. 8  Some  features  of  the  two  codes  of  interest  in 
the  current  study  are  described  below. 


^Dash,  S.M.,  and  H.S.  Pergament,  "A  Computational  Model  for  the 
Prediction  of  Jet  Entrainment  in  the  Vicinity  of  Nozzle  Boattails  (The 
BOAT  Code),"  NASA  Langley  Research  Center,  1978.  (NASA  CR-3175, 
publication,  UNCLASSIFIED.) 

8Fishburne,  E.S.,  and  A.K.  Varma,  "Investigations  of  Chemical  Reactions 
in  a  Turbulent  Media,"  Acta  Astronautics,  Vol.  6,  1979,  pp.  297-308. 
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BOAT  COMPUTER  PROGRAM 


The  BOAT  code  solves  a  parabolic  axisymmetric,  steady  viscous  flow 
in  which  the  turbulence  i3  described  by  a  two-equation,  turbulence 
kinetic  energy  model.  The  code  is  capable  of  handling  multi-species, 
multi-step  chemistry.  However,  it  does  not  calculate  the  correlations 
between  the  various  fluctuating  scalar  quantities  such  as,  species 
concentrations  and  between  species  concentrations  and  temperature  that 
are  expected  to  be  important  in  turbulence-chemistry  and 
turbulence-radiation  interactions.  The  program  calculates  only  two 

turbulence  _ quantities: _  the  turbulent  kinetic  energy, 

1/2q2  =  1/2(u' 2  +  v*  2  +  w" 2) ,  and  the  dissipation  rate,  e.  The 
chemistry  calculations  in  BOAT  neglect  the  effect  of  the  turbulent 
fluctuations  on  the  reactions,  and  only  the  mean  species  and 
temperatures  are  used  to  calculate  the  chemical  source  terms. 


RSL  COMPUTER  PROGRAM 

The  RSL  code  solves  a  parabolic,  two-dimensional  (planar  or 
axisymmetric) ,  steady,  viscous  flow  that  includes  the  conservation 
equations  for  the  mean  properties  and  all  the  independent  second-order 
correlations  for  a  simple  reacting  system.  The  code  has  been  designed 
to  handle  a  total  of  three  chemical  species  undergoing  a  one-step 
forward  reaction,  a  +  (3  y*  The  use  of  the  basic  shear  layer 
assumptions  leads  to  a  set  of  23  independent,  coupled, _  partial 
differential  equations;  equations  for  the  mean  variables,  p,  u^,  H,  a, 
and  p,  and  the  second-order  correlations.  The  computer  program  actually 
solves  a  total  of  29  equations  for  a  three- species  system.  The 
additional  equations  provide  a  check  on  the  mass  conservation  in  the 
program  and  on  the  accuracy  of  the  numerical  scheme. 

A  major  objective  of  the  RSL  code  development  was  to  provide  a 
research  tool  to  study  the  interactive  effects  of  turbulence  and 
chemistry.  In  the  current  application,  the  effect  of  all  the 
second-order  correlations  on  the  chemical  source  term  for  the  one-step 
reaction  can  be  included  in  the  calculation.  The  interaction  effects 
can  also  be  ignored  and  only  tte  mean  variables  used  in  the  chemistry 
calculation.  A  comparison  of  the  two  procedures  will  provide  a  measure 
of  the  importance  of  the  turbulence-chemistry  interaction  for  the  flow 
under  study. 


COMPARISON  OF  COMPUTER  PROGRAMS 

The  two  programs  have  different  strengths  and  weaknesses  as  far  as 
the  objectives  of  the  current  study  are  concerned.  The  objective  is  to 
evaluate  the  effect  of  turbulent  fluctuations  of  species  and  temperature 
on  the  mean  gaseous  radiation  from  an  afterburning  plume.  The  BOAT  code 
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can  handle  the  multi-step  chemistry  of  the  large  number  of  species  in  a 
tactical  missile  plume  with  the  neglect  of  turbulence-chemistry 
interaction  effects.  However,  the  BOAT  code  is  unable  to  directly 
provide  the  values  of  various  scalar  correlations  of  species  and 
temperature  in  the  flow  which  have  to  be  known  to  calculate  their  effect 
on  the  radiation.  In  the  two-equation  turbulence  model  approach,  these 
scalar  correlations  will  have  to  be  estimated  in  some  consistent  manner 
from  1/2q2  and  e.  The  RSL  code  solves  modeled  equations  for  these 
scalar  correlations  and,  therefore,  directly  provides  their  values 
across  the  flowfield,  without  the  need  for  additional  assumptions. 
However,  the  RSL  code  has  been  designed  for  a  simple  one-step  chemical 
reaction.  A  modeling  approach  for  multi-species,  multi-step  chemistry 
has  been  formulated  for  use  within  the  framework  of  the  RSL  code,2  but 
the  approach  requires  many  assumptions,  and  has  not  been  extensively 
tested.  It  was  decided  that  the  salient  features  of 
turbulence-chemistry  and  turbulence-radiation  interactions  can  be 
studied  in  a  simple  chemically  reacting  flow  that  models  one  aspect  of 
afterburning  in  a  real  plume.  The  reaction  parameters  and  flow 
conditions  were  chosen  to  simulate  properties  of  the  actual  plume  of 
interest.  With  the  use  of  a  simple  one-step  chemical  reaction,  it  is 
possible  to  compare  the  results  of  the  BOAT  and  RSL  codes  for  the  same 
flow  conditions. 


COMPARISON  OF  RESULTS  OF  THE  TWO  CODES  FOR  NONREACTING  FLOWS 

We  emphasize  that  the  comparison  of  predictions  from  these  two  codes 
was  not  a  primary  objective  of  this  study.  The  decision  wa3  made  to  use 
the  RSL  code  to  provide  both  the  mean  and  fluctuating  properties  of  the 
plume  because  we  believe  that  to  be  the  best  available  code  that 
provides  all  these  quantities.  In  addition,  it  was  an  initial  objective 
to  provide  predictions  from  the  current  codes  in  use  throughout  the 
industry  for  comparison,  for  both  the  flowfield  and  radiation.  An 
initial  calculation  of  a  turbulent  reacting  axisymmetric  jet  of  carbon 
iiionoxide  (CO)  in  a  coflowing  air  stream  with  the  BOAT  and  RSL  codes 
showed  very  significant  differences  between  the  results  from  the  two 
codes.  The  parameter  of  the  flow  and  the  results  are  discussed  later  in 
this  section.  The  differences  between  the  codes  for  the  reacting  flow 
necessitated  a  more  extensive  comparison  of  the  two  codes  for  simpler 
flows,  and  a  nianber  of  such  flowfield s  were  studied  in  order  to  isolate 
the  reasons  for  these  differences  and  also  to  compare  the  predictions 
with  measurements.  The  results  of  these  tests  are  discussed  below.  We 
include  this  comparison  here  because  it  serves  at  least  two  useful 
purposes.  It  establishes  the  validity  of  the  mean  flow  predictions  from 
RSL  for  this  type  of  plume  by  comparison  with  the  results  from  the  more 
familiar  SPF/1.  Secondly,  the  overall  validity  of  both  codes  is 
assessed  by  comparison  with  experiment. 
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Constant  Density  Axisymmetric  Jet  Into  Still  Air 

The  data  selected  for  this  test  was  that  obtained  by  Maestrello  and 
McDaid,  and  was  Test  Case  6  of  the  NASA  Free  Turbulent  Shear  Flows 
Conference.9  The  initial  mean  velocity  profile  was  measured  at  a 
location  1  jet  diameter  (x/r^  =  2)  downstream  of  the  exit,  and  was  U3ed 
to  start  both  the  codes.  The  BOAT  code  used  the  standard  k-e* 
turbulence  model.  The  q2  and  e  profiles  at  the  initial  station  are 
derived  with  the  use  of  a  mixing  length  model.  The  RSL  code  is  normally 
started  with  a  cosine  shaped  profile  for  the  turbulence  kinetic  energy, 
q2/2,  with  an  arbitrary  maximum  value  of  .003(uj  -  uj2.  The  length 
scale,  A.  is  related  to  the  shape  of  the  q2  profile;  for  an  axisymmetric 
jet  it  is  calculated  as  a  constant  times  the  distance  from  the  jet  axis 
to  the  point  where  q2  is  25%  of  its  maximum  value.  The  standard  value 
of  the  constant  is  0.2.  The  results  from  the  two  codes  are  shown  in 
Figure  1  along  with  the  experimental  results  for  the  decay  of  the  jet 
velocity  on  the  centerline.  Predictions  from  both  the  codes  are 
reasonably  satisfactory  for  this  flow.  The  RSL  code  predicts  a  shorter 
core,  but  shows  better  agreement  farther  downstream. 

The  results  for  the  two  codes  shown  in  Figure  1  are  based  on 
different  initial  turbulence  levels  and  scales.  It  is  feasible  to  match 
the  initial  turbulence  levels  and  profiles  in  the  two  cedes,  but  in  the 
current  formulation  it  is  not  possible  to  match  the  turbulence  scales 
over  the  entire  flowfield.  Currently,  RSL  does  not  incorporate  a 
transport  equation  for  the  scale.  The  scale  is  algebraically  related  to 
the  local  q2  profile,  and  is  constant  across  the  profile  at  every  axial 
station.  A  number  of  RSL  runs  were  made  after  matching  the  mean 
velocity  and  turbulence  kinetic  energy  profiles  to  the  results  obtained 
from  the  BOAT  code  at  x  =  Urj.  The  e  profile  was  used  to  estimate  the 
average  turbulence  scale  [ (q2/2) 3/2/e]  and  indicated  the  proper  constant 
to  be  ~0.13  at  the  initial  station.  It  must  be  remembered  that  the 
value  of  the  scale  constant  will  have  to  be  changed  at  other  stations  to 
continue  to  match  the  scale  being  used  in  the  BOAT  code,  and  this  cannot 
be  accomplished  in  the  current  RSL  code.  The  RSL  runs  used  three  scale 
constants,  0.1,  0.15,  and  0,2.  The  results  are  shown  in  Figure  2.  The 
results  show  that  by  selecting  a  smaller  scale  constant  that  is  in  line 
with  the  scale  in  the  BOAT  code,  the  RSL  code  results  .or  the  length  of 
the  core  are  in  good  agreement  with  the  BOAT  code.  Farther  downstream, 
it  appears  that  the  scale  constant  should  increase  towards  0.2,  the 
nominal  value  used  in  standard  RSL  runs. 


®Free  Turbulent  Shear  Flows,  Vols.  I  4  II,  NASA  Langley  Research  Center, 
July  1972.  (NASA  SP-321,  publication  UNCLASSIFIED.) 

*The  common  terminology  for  the  turbulence  kinetic  energy  is  k  =  1/2q2. 
Two-equation  turbulence  model  codes  are  therefore  referred  to  as  k-e 
models. 
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.1 —  .  BOAT,  k-c  turbulence  model 
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FIGURE  1.  Predictions  and  Measurements  of  the  Axis  Velocity 
in  a  Constant  Density,  Low  Speed,  Inert  Jet  Exhausting  into 
Still  Air.  NASA  Test  Case  6,  Uj  =  211  m/sec,  r^  =  .031  m. 
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FIGURE  2.  Effect  of  Varying  the  RSL  Turbulence  Length 
Scale  Constant  on  Axial  Velocity  Decay.  Initial  Turbulence 
Level  Matched  to  BOAT  Run  at  x  =  4rj.  Constant  Density, 

Low  Speed, Inert  Jet  Exhausting  into  Still  Air. 
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In  summary,  both  the  codes  are  in  reasonably  good  agreement  with  the 
experimental  data  for  low  speed  isothermal  jets  when  they  each  use  their 
normal  start  procedures  for  the  turbulence  parameters.  The  results  of 
the  two  codes  also  appear  to  be  consistent  with  one  another  when  the 
basic  differences  in  the  length  scales  used  in  the  programs  are  taken 
into  account. 


Heated  Axisymmetric  Jet  Into  Still  Air 

Measurements  of  the  decay  of  the  axial  velocity  in  heated 
axisymmetric  jets  into  still  air  are  reported  in  Abramovich. 10  The 
initial  ratio  of  jet  temperature  to  the  ambient  is  3.25.  Both  the  code3 
were  started  with  basically  top  hat  mean  profiles  at  the  exit  plane  of 
the  jet.*  The  length  scale  A  is  the  same  as  in  the  isothermal  jet 
study;  0.2  times  the  distance  from  the  jet  axis  to  the  point  where  q2 
drops  to  25%  of  its  maximum  value.  This  procedure  will  overestimate  the 
scale  to  some  degree  in  the  potential  core  region  where  the  turbulence 
is  restricted  to  the  shear  layer  region,  but  the  scale  specification 
will  be  correct  farther  downstream.  The  results  from  the  two  codes  are 
shown  in  Figure  3.  The  RSL  results  are  in  better  agreement  with  the 
data  for  this  flow.  The  BOAT  code  predicts  a  much  smaller  effect  of  the 
lower  density  of  the  jet  fluid  compared  to  the  results  for  the 
isothermal,  constant  density  flow. 

Th}  major  reason  for  the  differences  in  the  results  obtained  from 
the  two  codes  for  this  flow  is  the  different  choices  for  the  initial 
turbulence  intensity  and  length  scale.  Some  of  the  differences  are  also 
due  to  different  modeling  of  density  fluctuation  terms  in  the  full 
second-order  closure  turbulence  model  and  the  two-equation  model.  Some 
RSL  runs  were  started  with  the  initial  conditions  for  the  mean  velocity, 
mean  temperature,  and  turbulence  kinetic  energy  matched  to  the  results 
from  the  BOAT  code  at  x  =  15.32rj,  after  the  end  of  the  potential  core. 
The  average  A  across  the  profile  at  this  axial  position  in  the  BOAT 
results  corresponds  to  the  use  of  a  scale  constant  of  0.25  in  RSL 
instead  of  the  normal  value  of  0.2  used  in  RSL  for  axisymmetric  jet 
problems.  RSL  runs  were  made  with  the  scale  constant  set  equal  to  0.2 
and  0.25,  and  the  results  are  shown  in  Figure  4.  The  results  of  these 
RSL  runs  with  initial  conditions  matched  to  BOAT  at  x  s  15.32r^  show 
much  better  agreement  with  the  downstream  results  obtained  from  BOAT, 
than  the  results  shown  earlier  in  Figure  3.  The  RSL  run  with  the  scale 
constant  of  0.25,  that  approximately  matches  the  initial  scales  in  the 
two  codes,  deviates  somewhat  from  the  BOAT  results,  while  much  better 

1 °Abr amov ich ,  G . N . ,  The  Theory  of  Turbulent  Jets,  The  MIT  Press, 
Cambridge,  MA,  1963,  Chapter  7. 

•The  initial  q2  profile  in  RSL  was  chosen  to  be  sharply  peaked  near  the 
jet  radius,  with  a  maximum  value  of  ,003(uj  -  uj2. 
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FIGURE  4.  Effect  of  Varying  Turbulence  Length  Scale  Constant 
in  a  Heated,  Low  Speed  Jet.  BOAT,  RSL  Matched  at  x/r*  *  15.32. 

T./T  =  3.25,  u.  =  7  m/sec. 
j  00  j 
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agreement  is  shown  by  the  run  with  scale  constant  of  0.2.  The 
differences  between  the  RSL  and  BOAT  runs  with  the  matched  initial 
scales  may  be  due  to  a)  differences  in  turbulence  models  in  the  codes, 

b)  lack  of  agreement  among  the  scales  at  downstream  locations,  cr 

c)  differences  due  to  scale  variation  across  the  profile  in  BOAT,  while 
RSL  currently  uses  a  constant  scale  across  the  profile  for  a  jet  flow. 

In  summary,  the  RSL  code  demonstrates  significantly  better  agreement 
with  the  experimental  data  when  it  uses  its  normal  start  procedure.  The 
BOAT  code  shows  a  very  small  effect  of  the  density  variation  on  the 
axial  velocity  decay  profile,  and  is  not  in  good  agreement  with  the 
data. 


Hydrogen  Axisynanetric  Jet  Into  A  Moving  Air  Stream 

Eggers9  made  measurements  of  mixing  in  relatively  high  speed  streams 
of  hydrogen  and  air.  The  centerline  variation  of  the  axial  velocity  and 
hydrogen  mass  fraction  was  measured.  The  initial  velocity  profile  at 
the  nozzle  exit  shows  the  presence  of  wall  boundary  layers,  and  has  a 
substantial  velocity  defect.  The  RSL  code  has  some  difficulties 
starting  from  this  initial  profile  due  to  the  problem  of  specifying  the 
proper  turbulence  scale  distribution.  A  BOAT  run  was  carried  out  from 
the  initial  profile,  and  then  an  RSL  run  was  started  using  the  results 
of  the  BOAT  code  at  an  axial  station  4  radii  downstream  from  the  nozzle 
exit.  The  RSL  run  matched  the  mean  velocity,  mean  species 
concentration,  and  the  turbulent  kinetic  energy  profile  at  the  initial 
station.  The  average  turbulent  length  scale  across  the  profile  at  this 
axial  station  in  the  BOAT  results  corresponds  to  the  use  of  a  scale 
constant  of  0.2  in  RSL.  The  axis  velocity  and  hydrogen  mass  fractions 
obtained  from  the  BOAT  and  RSL  codes  (C  =  0.2)  are  shown  in  Figures  5 
and  6,  respectively  along  with  the  experimental  data.  The  BOAT  code 
predictions  are  in  good  agreement  with  the  data.  The  RSL  results  start 
out  in  good  agreement  with  the  BOAT  results  (as  expected  for  matched 
initial  conditions)  but  then  deviate  from  them.  We  compared  the  length 
scales  in  the  two  codes  at  x/rj  =  40  and  found  that  the  average  scale 
across  the  flow  in  BOAT  was  50%  higher  than  that  in  RSL.  Farther 
downstream,  at  x/rj  =  100,  the  scales  in  the  two  codes  are  again 
comparable.  We  made  another  RSL  run  using  a  scale  constant  C  =  0.3,  and 
these  results  a'  e  also  plotted  in  the  figures.  The  velocity  decay  data 
is  bounded  by  the  R£f.  results  with  the  two  length  scales.  The  species 
mass  fraction  is  over predicted  at  large  distances,  even  with  the  larger 
scale  constant.  It  is  clear  from  these  comparisons  that  much  of  the 
difference  between  the  mean  flow  predictions  of  these  two  codes  for  this 
particular  test  case  is  due  to  the  different  length  scales.  Due  to 
limitations  in  the  length  scale  model  currently  being  used  in  RSL,  the 
RSL  results  for  this  flow  are  not  very  satisfactory.  The  results  of  tne 
BOAT  code  for  this  flow  are  in  good  agreement  with  the  data. 
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FIGURE  5.  Axis  Velocity  in  a  Hydrogen  Jet  in  a  Coflowing 
Stream.  NASA  Test  Case  12,  u\  =  1074  m/sec,  r.  =  .0058  m, 
uoo  =  394  m/sec.  3  J 
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Summary  of  RSL  -BOAT  Compari sons  For  Non-Reacting  Flows 

~  The  three  test  cases  carried  out  here  show  mixed  results.  Both  the 
codes  are  satisfactory  for  the  constant  density  jet  flow,  but  for  the 
two  variable  density  flows,  RSL  does  better  for  the  heated  jet 
experiments  while  BOAT  performs  better  for  the  light  species  jet. 
Further  tests  on  more  complex  flows  will  have  to  be  carried  out  to 
document  the  superiority  of  one  code  over  the  other.  The  predictions  of 
the  codes  are  in  reasonable  agreement  with  each  other  when  the  initial 
conditions  for  the  turbulence  scale  and  turbulent  kinetic  energy  are 
matched . 


REACTING  FLOW  CALCULATIONS 

The  two  computer  codes  were  used  to  calculate  the  flowfield 
properties  of  a  turbulent  reacting  axisymmetric  jet  of  carbon  monoxide 
in  a  coflowing  stream  of  air.*  The  initial  conditions  for  the 
calculations  were  as  follows: 


TABLE  1.  Initial  Conditions  for  Model 
Afterburning  Plume 


Jet  Stream 

. 

External  Stream 

Uj  s  2500  rn/  sec 

u*  =  300  m/sec 

D  =  0.13  m 

T.  =  300°K 

Tj  =  1200°K 

Species  Mass  Fraction 

Species  Mass  Fraction 

CO  =  0.25 

02  =  0.23 

N2  =  0.75 

N2  =  0.77 

^chosen  to  simulate  the  exhaust  plume  from  a  small  <5000  lbf  thrust) 
tactical  missile  using  a  non- alum ini zed  solid  propellant. 
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For  both  oases,  the  initial  velocity  a'd  species  profiles  were 
cosine  curves  between  the  jet  and  external  stream  values.  The  initial 
temperature  profile  had  a  peak  temperature  of  1836  K  in  the  mixing 
region  to  ignite  the  CO  jet.  The  BOAT  code  used  its  conventional  start 
procedure  (a  mixing  length  model)  to  calculate  the  initial  profiles  of 
turbulence  kinetic  energy,  1/2q2,  and  the  dissipation  rate,  e.  The  BOAT 
code  incorporates  an  empirical  compressibility  correction  model  to  be 
used  for  high  speed  flows,  and  the  calculations  for  the  CO  jet  were  made 
both  with  and  without  the  use  of  this  correction.  The  RSL  code  was  also 
started  with  its  normal  start  procedure.  A  cosine  shaped  profile  for 
the  turbulence  kinetic  energy  1/2q2  is  used  with  a  maximum  value  of 
,003(uj  -  u*)2.  The  length  scale,  A  is  based  on  the  q2  profile  with  a 
scale  constant  of  0.2. 

The  carbon  monoxide  reaction  with  atmospheric  oxygen  is  assumed  to 
be  an  overall  one-step  irreversible  reaction 


a  + 


Y 


(3) 


with  a  the  air  species  (N2,02), 
p  the  plume  species 
Y  the  products  (C02,N2). 


k  =  2  x  1013  exp(-15000/RT) 


(4) 


This  reaction  mechanism  is  a  significant  simplification  of  the 
multi-step  elementary  reactions  that  are  believed  to  take  place.  The 
chemical  source  terms  in  the  BOAT  code  are  calculated  using  the  mean 
species  concentrations  and  the  mean  temperature,  and  neglecting  the 
effect  of  turbulent  fluctuations  on  the  chemistry.  RSL  calculations 
were  carried  out  in  two  different  ways;  one  using  only  the  mean  values, 
and  the  other  including  the  effect  of  second-order  species  fluctuation 
correlations  and  mean  square  temperature  fluctuations  on  the  chemical 
source  terms. 


The  axial  variation  of  the  maximum  temperature  in  the  flow  obtained 
from  the  BOAT  and  RSL  computer  programs  is  plotted  in  Figure  7.  The  use 
of  the  compressibility  correction  in  the  BOAT  predictions  decreases  the 
rate  of  mixing  and  extends  the  reaction  zone  to  some  degree.  The  RSL 
calculations  shown  here  used  only  the  mean  scalar  variables  to  calculate 
the  chemical  source  terms,  that  is,  the  same  procedure  as  in  the  BOAT 
code.  The  predictions  of  the  RSL  r.»n  with  the  normal  value  of  the  scale 
constant,  C  =  0.2  for  axisymmetric  jet  flows,  are  significantly 


22 


NWC  TP  6307 


FIGURE  7.  Peak  Temperature  in  the  Model  After¬ 
burning  Exhaust  Plume.  Temperature  Normalized  by 
a  Reference  Temperature  of  300  K. 
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different  from  both  the  BOAT  predictions.  A  comparison  of  the  values  of 
the  turbulence  scale  in  the  two  codes  at  different  axial  stations  was 
carried  out,  and  the  results  are  tabulated  below. 


TABLE  2  -  Average  Turbulence  Scale  Constants,  A/2r j 


BOAT 

RSL 

scale  constants 

x/rj 

with  corap 

w/o  corap 

C=0. 15 

mm 

C=0.25 

0 

mm 

0.10 

mm 

0.17 

0.21 

20 

0.37 

wmmmm 

0.27 

40 

■a 

0.58 

0.19 

■ 

0.47 

100 

1.08 

0.87 

0.36 

0.72 

1.12 

At  the  initial  station,  the  RSL  scale  length  is  larger,  but  at  other 
axial  stations,  within  the  afterburning  region  (x/rj  <  100),  the  scale 
length  is  significantly  smaller.  At  around  x/rj  =  100,  the  difference 
becomes  less  than  at  earlier  stations.  Additional  RSL  runs  were  made 
with  the  scale  constant  set  at  0.15  and  0.25  to  attempt  to  resolve  the 
discrepancy  between  the  RSL  and  BOAT  results.  The  complete  set  of 
results  for  the  RSL  cases  is  shown  in  Figure  8. 

The  scale  in  RSL  runs  with  C  =  0.25  is  closer  to  the  BOAT  scale  in 
the  region  from  10  to  50  diameters,  and  so  are  the  predictions  for  the 
peak  temperature.  Still,  there  are  other  differences  between  the  two 
results  that  cannot  be  explained  by  the  different  scales,  but  certainly 
matching  the  scales  improves  the  agreement  between  the  predicted  results 
from  the  two  codes. 

For  the  purposes  of  the  current  study,  we  need  the  consistent  set  of 
species  and  temperature  correlations  provided  by  the  RSL  code.  In  the 
absence  of  experimental  data  on  reacting  afterburning  plume  flows  to  fix 
the  empirical  value  of  the  scale  constant,  we  will  retain  the  standard 
value  of  0.2  that  we  use  for  nonreacting  axisymmetric  jet  flows. 


TURBULENCE-CHEMISTRY  INTERACTION 

The  non-linear  dependence  of  the  chemical  reaction  rate  on 
temperature  [Eq.  (4)3  suggests  that,  like  the  radiant  emission,  the 
time-averaged  chemical  species  production  is  different  from  that 
evaluated  at  the  mean  temperature  and  species  concentrations. 
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FIGURE  8.  Peak  Temperature  in  the  Model  After¬ 
burning  Exhaust  Plume  for  Various  Length  Scale 
Constants.  Peak  Temperature  Normalized  by  a 
Reference  Temperature  of  300  K. 
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Expansions  of  the  species  production  rate,  keeping  terms  up  to  the 
second-order  correlations  are  included  in  RSL.  The  rate  of  formulation 
of  product  species,  y, .is 


dy/dt  =  k|>  a  6 


(5) 


If  each  term  is  expanded  in_terms  of  a  mean  and  fluctuating  part, 

kf  =  kf  +  kj>,  a  =  a  +  ct1  ,  p  =  p  +  p'  ,  then 

dy/dt  =  kfa  g(  1  +  a'p'/ap+  . .  . )  .  (6) 

where 

kf  =  kf(T)[  1  +  E(E/2  -  1)F2/T2  +...].  (7) 


The  neglected  terms  involve  second-order  correlations  of  kf  with  a' ,  p' , 
and  all  third-  and  higher-order  correlations.  a'  p'/a  p  is  the 
"unmixedness  correlation”  which  defines  the  effect  of  species 
fluctuations  on  the  reaction  rate.  In  *  diffusion  flame  where  the 
reactants  are  initially  separated  (as  the  case  for  afterburning  rocket 
plume),  this  correlation  is  negative.  Limiting  cases  are:8 


fast  chemistry:  a'p'  /a  p  +  -1.0 


slow  chemistry:  a'p'  /j  p  +  -0.2  . 


RSL  runs  have  been  made  using  the  above  expression  for  the  chemical 
source  term  _asi  well  as  using  only  the  mean  variables;  that  is, 
dy/dt  =  kf(T)a  p.  The  two  results  are  shown  in  Figure  9.  The 
calculations  show  that  in  this  simulation  of  an  afterburning  plume,  the 
chemical  reactions  are  virtually  completed  at  50  radii  downstream  of  the 
nozzle  exit.  The  fuel  (CO)  mass  concentration  has  decreased  from  an 
initial  value  of  0.25  to  about  0.02  at  this  station.  Further 
downstream,  the  flow  is  simply  the  mixing  and  cooling  of  a  hot 
axisymmetric  jet.  The  turbulence-chemistry  interaction  effects  will 
only  be  important  in  the  initial  reacting  region,  and  the  inclusion  of 
the  species  fluctuations,  as  expected,  acts  to  reduce  the  reaction  rate, 
and  leads  to  a  slightly  lower  peak  temperature.  Due  to  the  small  extent 
of  the  reaction  zone,  turbulence-chemi3try  interactions  are  not  very 
important  in  this  flowfield.  The  axial  profile  of  the  peak  temperature 
is  virtually  identical  for  the  two  runs  beyond  50  radii. 
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FIGURE  9.  Effect  of  Turbulent  Fluctuations  on 
Peak  Temperature  in  the  Model  Afterburning 
Exhaust  Plume.  Peak  Temperature  Normalized 
by  a  Reference  Temperature  of  300  K. 
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The  fact  that  the  turbulence-chemistry  interaction  is  not  a  major 
factor  in  this  particular  flow  can  be  explained  by  two  effects.  One  is 
that  the  magnitude  of  the  "unmixedness  correlation"  is  insufficient  to 
substantially  reduce  the  reactant  production  below  its  nominal 
time-averaged  value.  Moreover,  the  effect  of  the  temperature 
fluctuations  acts  in  the  opposite  direction  to  increase  the  rate 
constant  (for  exothermic  reactions).  These  two  effects,  unmixedness  and 
temperature  fluctuations,  therefore  tend  to  cancel  each  other.  It  is 
important  to  note  also  that  this  result  is  highly  dependent  on  the 
specific  reactions  because  the  augmentation  by  temperature  fluctuations 
$epends_  on  activation  energy.  For  the  case  considered  here, 
E  =  E/RT  <7.5  in  the  reacting  layer  and  T'2/T2  <0.1.  For  these 
conditions,  the  rate  constant  is  at  most  tripled,  which  is  insufficient 
to  significantly  change  the  mean  flow  properties  even  if  there  were  no 
compensating  unmixedness  correction.  The  detailed  results  for  the 
exhaust  plume  certainly  support  these  qualitative  arguments.  We  caution 
that  each  flow  situation  and  chemical  system  must  be  individually 
examined  for  the  relative  importance  of  these  effects.  More 
importantly,  the  conclusions  here  are  based  on  a  global  reaction  scheme. 
It  is  possible  that  a  more  precise  specification  of  the  chemistry,  in 
terms  of  the  elementary  reactions  involving  all  the  trace  species,  will 
be  more  sensitive  to  the  turbulence. 


COMPARISON  OF  TURBULENCE  QUANTITIES 

The  final  comparison  between  the  RSL  and  BOAT  code  predictions  is 
between  the  radial  profiles  of  mean  temperature,  and  the  turbulence 
quantity  predicted  by  both  codes;  the  turbulence  kinetic  energy.  The 
mean  square  temperature  fluctuation  T'  2  and  product  species  fluctuation 
yrz*  from  the  RSL  code  are  also  shown.  The  radial  mean  temperature 
profiles  at  x/rj  =  40  and  x/rj  =  100  and  radial  profiles  of  the 
turbulence  kinetic  energy,  temperature  and  species  correlations  are 
plotted  in  Figures  10  and  11.  The  profiles  for  the  mean  temperature  and 
the  turbulence  kinetic  energy  at  the  same  axial  stations  from  the  BOAT 
code  are  3howr»  in  Figures  12  and  13  without  and  with  the  compressibility 
correction,  respectively.  Due  to  the  extended  reaction  zone  in  the  BOAT 
predictions  with  the  compressibility  correction,  the  mean  temperature 
profiles  at  x/rj  =  40  are  quite  different.  The  turbulence  level  with 
the  compressibility  correction  is  increased  over  that  without  the 


•This  quantity  is  a  mean  square  mass  fraction  fluctuation.  The 
radiation  formulation  requires  mole  fractions.  The  conversion  from  mass 
fraction  to  mole  fraction  in  this  case  is  c*  =  y'  c/y,  approximately. 

In  the  most  general  case,  conversion  from  fluctuating  mass  fraction  to 
mole  fraction  involves  correlations  between  all  species.  The  largest 
error  in  neglecting  the  fluctuations  occurs  in  mixtures  with  species 
with  disparate  molecular  weights  (e.g.,  hydrogen  containing  mixtures). 
The  error  is  negligible  for  the  mixture  considered  here. 
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FIGURE  10.  Radial  Distribution  of  Mean  and  Flue 
tuating  Properties  Predicted  by  RSL  Code  for  the 
Model  Afterburning  Exhaust  Plume  at  x/rj  *  40. 
Temperature  and  Velocity  Normalized  by  Reference 
Values  of  300  K  and  2500  m/sec,  Respectively. 


FIGURE  11.  Radial  Distribution  of  Mean  and  Fluctu¬ 
ating  Properties  Predicted  by  RSL  Code  for  the  Model 
Afterburning  Exhaust  Plume  at  x/rj  =  100.  Temperature 
and  Velocity  Normalized  by  Reference  Values  of  300  K 
and  2500  ra/sec.  Respectively. 
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FIGURE  12.  Radial  Profiles  of  Mean  and  Fluctu¬ 
ating  Properties  Predicted  by  the  BOAT  Code 
without  a  Compressibility  Correction  for  the 
Model  Afterburning  Exhaust  Plume  at  x/rj  =  40 
and  100.  Temperature  and  Velocity  Normalized  by 
Reference  Values  of  300  K  and  2500  m/sec. 
Respectively. 


31 


2-  T  / 

X=  lOOrj 


\  -1.2 

\ 

\ 


0  I  2  3.4  5  6 

y/rj 

FIGURE  13.  Radial  Profile  of  Mean  and  Fluctu¬ 
ating  Properties  Predicted  by  the  BOAT  Code  with 
a  Compressibility  Correction  for  the  Model  After¬ 
burning  Exhaust  Plume  at  x/rj  =  40  and  100. 
Temperature  and  Velocity  Normalized  by  Reference 
Values  of  300  K  and  2500  m/sec,  Respectively, 
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correction,  and  is  larger  than  that  predicted  by  the  RSL  code. 

'  The  turbulent  kinetic  energy  levels  predicted  by  the  RSL  code  are 
comparable  to  the  levels  in  the  BOAT  predictions  without  the 
compressibility  corrections  in  spite  of  the  differences  in  the  mean 
properties.  This  agreement  i3  not  a  verification  of  the  predictions 
from  either  code,  but  it  is  an  encouraging  result,  particularly  since 
the  other  turbulence  correlations  predlutc'1  by  RSL  are  to  be  used  to 
study  the  effect  of  fluctuations  on  the  radiatio*'  from  the  hot  reaction 
products  without  benefit  of  a  similar  comparison  The  differences  in 
the  predictions  of  the  mean  properties  in  the  two  codes  for  reacting 
flows  require  further  investigation  and  detailed  comparison  with 
experimental  measurements  to  assess  the  relative  performance  of  the  two 
codes. 

The  RSL  results  for  the  scalar  correlations,  y' 2  and  T ' 2 ,  and  the 
velocity  correlation  a2,  show  that  these  radial  profiles  are  quite 
different.  It  is  not  correct  to  directly  relate  these  other  scalar 
correlations  to  the  turbulent  kinetic  energy  alone.  Thus,  the  coupling 
of  turbulence  and  chemistry  and/or  turbulence  and  radiation  that 
requires  information  about  many  second-order  correlations,  cannot  be 
adequately  studied  with  a  flowfield  code  using  simply  the  two-equation, 
k-e  turbulence  model  alone.  Additional  turbulence  equations  and  model 
assumptions  will  have  to  be  used  for  this  task.  For  example, 
three-equation,  k-e-g  turbulence  models,  algebraic  3tress  modeling 
(quasi-equilibrium  models) ,  or  complete  second-order  closure  models  as 
used  in  the  RSL  code  are  required. 


RADIATION  FROM  A  FLUCTUATING  FLAME 

We  begin  with  the  integrated  equation  of  transfer  for  the  spectral 
radiance  along  a  given  optical  path  through  the  plume 


kw  is  the  spectral  absorption  coefficent  which  is  due  to  the  summation 
of  the  contribution  from  many  rotational  lines  in  a  vibration-rotational 
band.  The  spectral  dependence  of  this  absorption  is  very  complicated; 
it  varies  by  an  order  of  magnitude  within  spectral  intervals  comparable 
to  a  line  width  (<  .05  cm"1),  and  there  may  be  hundreds  of  lines  of 
varying  strength  contributing  to  the  absorption  at  a  given  wavenumber. 
In  the  typical  spectral  resolution  of  low-resolution  radiance 
predictions  (5  -  25  cm"1),  there  may  be  contributions  from  thousands  of 
lines.  There  are  two  ways  of  dealing  with  this  computational 
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problem:  one  is  to  perforin  the  calculation  on  a  purely  monochromatic 
basis  (line-by-line  approach)  and  accept  the  computational  penalty;  the 
other  is  to  use-  a  band  model  approach  in  which  the  absorption  is 
spectrally  averaged  over  many  lines  before  the  spatial  integration  is 
performed.  For  homogeneous  optical  paths,  this  interchanging  of  the 
spectral  and  spatial  integration  is  exact.  However,  for  non-homogeneous 
paths,  it  is  exact  only  for  optically  thick  or  optically  thin  paths. 
Intermediate  optical  depths  are  essentially  an  inter pvlation  between 
these  two  limits  that  requires  spatial  averages  of  two  parameters:  the 
absorption  coefficient,  and  the  ratio  of  average  line  width  to 
spacing.11*12  These  spatial  averages  make  application  of  the  band  model 
calculation  difficult  to  apply  when  the  optical  path  is  to  include  the 
fluctuating  quantities.*13  The  line-by-line  approach  was  therefore 
chosen  for  this  study.  This  permits  the  direct  expansion  of  the 
radiance  in  terms  of  known  mean  and  second-order  correlations  of  the 
fluctuating  properties  provided  by  the  turbulent  flowfield  solution. 

There  is  a  disadvantage  to  choosing  the  line-by-line  formulation. 
The  required  tabulation  of  line  properties  (spectral  positions, 
strengths,  and  quantum  states)  is  incomplete.  Vibrational  transitions 
between  energy  states  that  beoome  populated  only  at  high  temperatures 
(>  600°k)  are  absent  from  the  available  compilation. 14  However,  an 
alternative  source  for  these  "hot  lines"  has  been  devised  for  the  linear 


1 1 Lud wi g ,  C.B. ,  et  al.f  Handbook  of  Infrared  Radiation  from  Combustion 
Gases,  1973*  (NASA  SP-3080,  publication  UNCLASSIFIED.) 

12y0ung,  S.J.,  "Nonisothermal  Band  Model  Theory,"  The  Aerospace 
Corporation,  September  1976.  (SAMS0-TR-76-207,  publication 

UNCLASSIFIED.) 

•Rhodes13  used  the  clever  approach  of  defining  an  ensemble  average  of 
the  radiance  over  many  realizations  of  the  optical  path.  Properties  at 
each  point  were  chosen  by  random  selection  from  an  assumed  probability 
density  distribution.  The  mean  radiance  was  obtained  by  arithmetrically 
averaging  the  radiance  for  each  realization  of  the  optical  path. 

13Rhodes,  R.P.,  "The  Effect  of  Turbulent  Fluctuations  on  the  Infrared 
Radiation  from  Rocket  Plumes,"  12th  JANNAF  Pliane  Technology  Meeting, 

CPIA  Publication  332,  Vol.  Ill,  December  1980,  pp.  81-97.  (paper, 
UNCLASSIFIED.) 

14McClatchey,  R.A.,  et  al.,  "AFCRL  Atmospheric  Absorption  Line 
Parameters  Compilation,"  Air  Force  Geophysics  Lab.,  January  1973. 
(AFCRL-TR-73-0096,  publication  UNCLASSIFIED.) 
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C02  molecule.15  In  the  optically  thin  limit,  the  band  model  and 
line-by-line  computations  are  equivalent.  The  temperature  sensitivity 
of  emission  from  the  4.3  COg  band  has  recently  been  described  in  this 
limit  for  the  band  model  formulation. 16  However,  the  need  in  the 
present  study  was  for  an  approach  that  was  uniformly  valid  for  any 
optical  depth.  Therefore,  the  line-by-line  formulation  was  chosen. 
Because  the  pressure  is  atmospheric  and  the  temperatures  are  moderate, 
we  assumed  a  Lorentz  line  shape  and  considered  only  collision 
broadening. 

A  precise  evaluation  of  the  radiance  from  a  fluctuating  optical  path 
requires  a  joint  probability  density  function,  PDF(T,  c^  St) ,  for  the 
temperature  and  composition  at  every  point  in  the  plume,  t,  along  the 
optical  path.  However,  the  properties  available  are  the  mean  values  T, 
cit  and  second-order  correlations,  T'^,  c'T' ,  and  cjei  at  each  point. 
The  fluctuating  radiance  calculated  here  can  use  only  these  quantities. 
The  approach  was  to  decompose  each  variable  into  its  mean  and 
fluctuating  part.  For  primary  variables,  these  two  components  are 
obtained  directly  from  the  turbulent  flowfield  solution.  Secondary 
quantities  such  as  the  Planck  function  and  absorption  coefficient  were 
expanded  about  their  values  at  the  mean  temperature  and  composition  to 
express  their  fluctuations  in  terms  of  the  fluctuating  species  and 
temperature. 

B  (T)  =  B  +  B'  s  B  (T)  +  T  dB/dT  T'/T  ♦  1/2T2  d2B/dT2  T'2/T2  +  ... 
worn  ^ 

k  (T)  =  k  +  k»  =  k  (T,c)  +  T  Bk/BT  T'/T  +  Bk/Bc  c' 

0J  o  o 

+  1/2  T2  a2k/BT2  T'2/T2  +  T  a2k/3T&c  c'T’/T  +  1/2  B2k/Bc2  c'2  ♦  ...  . 

(10) 


These  two  expansions,  along  with  T  =  T  +  T'  and  c  s  c  ♦  c'  were  used  in 
Eq.  8,  which  was  expanded,  retaining  only  terms  up  to  and  including  the 

t'  2  c ' T ' ,  and  c'2 


second-order  correlations  T 


The  algebra  is  tedious. 


l5McCullough,  R.W.,  "Compilation  of  CO2  Rotational  Lines  Including 
Transitions  That  Are  Important  at  Temperatures  up  to  2000°K, 
Aeronautical  Research  Associates  of  Princeton,  Inc.,  unpublished, 
December  1980.  (paper,  UNCLASSIFIED.) 

l^Draper,  J.S.,  L.S.  Bernstein,  and  W.K.  Cheng,  "Fluctuating  Emission 
from  Molecular  Vibrational-Rotational  Bands,"  J.  Quart.  Spectrosc.  and 
Radiat.  Transfer,  Vol.  23,  1980,  pp.  323-326. 
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and  only  the  final  result  is  given  here. 

The  particular  property  of  the  fluctuating  radiance  we  chose  to 
evaluate  was  the  time-averaged  value,  Nu(T,c).  This  quantity  is  what 
would  be  measured  by  a  sensor  that  ha3  a  temporal  response  that  is  slow 
compared  to  the  frequencies  of  the  plume  fluctuations.  There  are  other 
quantities  of  interest  which  define  the  effect  of  turbulent 
fluctuations,  such  as  the  mean  square  radiance  fluctuation,  (Nw  -  Nu)2. 
However,  this  quantity  does  not  indicate  the  direction  in  which  the 
fluctuations  aot  on  the  radiance;  enhancing  or  diminishing  the  average 
radiance  from  the  value  that  was  obtained  using  only  mean  flow 
properties. 

The  final  result  for  the  average  radiance  is  given  below.  The  steps 
taken  in  its  derivation  were: 


1)  expand  each  term  in  the  integral  in  terms  of  it3  mean  and 
fluctuating  component, 

2)  expand  the  transmittance  integral  about  the  mean  optical  path, 
separating  out  the  optical  depth  of  the  mean  flow;  and  then 
expand  the  fluctuating  part: 


(11) 


where  F(c’ ,T'...)  denotes  terms  involving  the  fluctuating  properties. 


Note  that  the  square  of  the  integral  within  the  expansion  can 
be  expressed  as  a  double  integral  over  the  square  domain  with 
area  s2. 

3)  Perform  the  multiplications  and  discard  all  higher-order  terms 
beyond  the  second-order  correlations, 

4)  time-average  the  result. 


In  the  time-average,  all  terras  involving  c'  and  T'  alone  have  zero 
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averages  (from  the  definition  of  the  fluctuations).  All  terms  involving 
products.  T’2,  o*T',  c* 2,  give  the  second-order  correlations  T'  c'T*, 
and  c'  2.  A  special  case  of  the  second-order  correlations  occurs  in 
evaluating  the  integral  terms.  The  following  forms  of  integrals  must  be 
averaged. 


A( s* )  T*  (s' )  B( s")  T'(s")  ds»  ds"  ds 


(12a) 


(12b) 


where  A(s),  B(s),  and  C(s)  are  terms  involving  derivatives  of  Bu  and  k 
evaluated  at  the  mean  flow  properties  at  the  point  s.  We  express  these 
time-averages  as. 


a 


■jCT 


s  I  us 

A(s)  T'(s)  B(s')  T'(s')  ds*  d3  =  |  |  A(s)  B(s')  T'(s)  T'(s')  ds*  ds 

(13a) 


f  C(s)  f  f  A(s’)  B(s" 

Jo  JoJo 


)  T*  (s' )  T'(s")  ds'  ds"  ds  = 


A  AA 

=  I  C(s)  II  A(s')  B(s" 

Jo  JcJo 


)  T'CsMT'ts")  ds1  ds"  ds 


(13b) 


These  averages  give  the  two-point  correlation  of  the  turbulent 
fluctuations  along  the  optical  path.  In  the  subsequent  analysis,  these 
correlations  are  expressed  in  terras  of  the  correlation  coefficients. 
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T’(s)  T* (s’ )  =  [T’2(s>] 


[T‘2(s»>3 


*rTT(3t3’) 


There  are  similar  correlations  between  c'  and  T* ,  and  between  c*  and  c’ 
In  all  cases,  we  use  an  exponential  fora  of  the  correlation 


rTT(s,s’)  =  exp(-  Is  -  s’ \/ij 


where  A  is  a  correlation  length  which  we  take  to  be  the  same  for  all  the 
scalar  correlations.  It  was  our  original  intent  to  choose  different 
functional  forms  for  the  correlation  to  test  the  sensitivity  of  the  mean 
radiance  to  these  terms.  However  ,  this  particular  form  admits  a 
convenient  transformation  that  gives  an  accurate  numerical  quadrature, 
and  so  this  form  was  maintained  throughout  the  study.  We  argue  that 
this  term  always  appears  as  an  integral  and  therefore,  the  sensitivity 
to  the  precise  form  of  the  integrand  is  reduced.  The  correlation 
length.  A,  is  a  parameter  that  was  varied  over  a  wide  range. 

Following  the  procedure  outlined  above  gives  the  equation  for  the 
time  averaged  spectral  radiance  from  a  turbulent  flow.  As  in  most 
radiative  transfer  problems,  we  have  taken  the  step  of  using  the  mean 
optical  path  length  rather  than  the  geometric  path  length  as  the 
variable  of  integration 


•jp 


i  ~  Ike  ds' 

L  o 


{-£  *°° ds) ds  =jtL 


F(s)  kQc  expl- I  kQc  ds’J  ds  =  /  F[s(t)]  e  dx  .  (17) 


The  full  equation  for  the  time-averaged  radiance  is 
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K  ■  jf L  v-'*  ♦ 

*  jfLfe  pWo  *  Vo  *  i  f‘Bir)fTI/f'  ’ 

+  ^TB0kT/ko  ♦  TBjk0/k0  ♦  TB0kT/koe  ♦  TBj/Z'j  c'T'/T 

♦  B0kcC/k0  ♦  B0kc/k0o^c,ije-TdT  ♦ 

jCL=«r  2  T^k-jo  T'^/T^  +^TkcTo  +  TkTjo'T’/T  ♦  ^kocc  ♦  kc^ 


(184) 


(18b) 


c'^  ds'  e~Tdt 


-J*L  (fBokT/k0  *  TB^T^/T^’^jf8  fkT0  ^T2/x^1/2  rTT(S( 


s’)  ds'  e-Tdt 


-  jTL  (tBoV^o  *  TBt^T^/T ♦  k0^^1/2  rcT(s,s')  ds'  e-^dx 

-  jfl  (w-o * v)H’,2£  TkTo^T,i/T^1/2  rcT(s,s')  ds'  e"Tdt 


-/L  (wo  •  vH'tf  (k=:  *  ‘.X"')’" 


r0C(3,s')  ds1  e_Tdx 


(18c) 


(18d) 


(18e) 


(18£) 


(18g) 


?/L  Bo  X%3[?‘<TC(T,Z/T^)1/2|S,’lTkTIf7T/T^  1/2U*rTT<s,»srt>ds'da"  ^ 


(18h) 


•t  JX  •  ‘.)Hl 


3,rcT(s' ,s")ds,ds"  e"*dx 


(181) 


r„_(s' ,s")ds,ds"  e“Tdx 
00 

(18  j) 


•We  use  the  shorthand  subscript  notation  for  derivatives,  for  example, 
Bt  =  dB/dT,  kTT  =  32k/3T2. 
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This  expansion  is  uniformly  valid  for  all  optical  depths,  but  it 
does  require  that  the  perturbations  in  the  optical  path  are  small 
compared  to  the  -mean  optical  path  because  of  the  expansion  of  the 
exponential  [Eq.  (11)].  The  error  in  this  expansion  vanishes  as  -*■  0 
and  as  <*>.  Numerical  tests  show  that  the  relative  error  is  maximum 
for  =  0(1),  but  that  it  is  at  most  2  percent  even  if  the 
perturbations  in  optical  depth  are  as  large  as  half  the  mean  optical 
depth.  The  perturbations  for  the  levels  of  fluctuations  examined  here 
were  much  smaller  than  that,  and  the  expansion  is  adequate  for  the 
present  application. 

The  physical  interpretation  of  the  terms  in  the  radiance  is  given 
below: 


Term  Physical  Description 

18(a)  Radiance  from  the  source  function  BQ,  evaluated 

at  the  mean  temperature,  and  transmitted 
through  the  optical  path  evaluated  at  the  mean 
properties. 

18(b)  Radiance  from  the  fluctuation  augmented  source 

function  transmitted  through  the  optical  path 
evaluated  at  the  mean  properties. 

18(c),  18(h)  -*•  18(j)  Radiance  from  source  function  evaluated  at  mean 

temperature  transmitted  by  fluctuating  optical 
path. 

18(d)  •>  18(g)  Radiance  from  fluctuating  source  function 

transmitted  by  fluctuating  optical  path. 


Another  viewpoint  that  may  clarify  the  expansion  is  to  view  each  term  as 
a  product  of  two  factors:  a  local  emission  (source)  term  and  a 
transmittance  term.  The  expansion  in  Eq.  18  can  then  be  schematically 
outlined  as  shown  below. 


Terms  in  Eq.  (18) 


Source 


Transmittance 


18(a)  mean 

18(b)  fluctuating 

18(c)  -►  18(h)  +  18(j)  +  18(k)  mean 

18(d)  +  18(e)  +  18(f)  +  18(g)  fluctuating 


mean 

mean 

fluctuating 

fluctuating 
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Each  of  these  integral  terms  was  evaluated  numerically  for  the  known 
distribution  of  mean  and  fluctuating  properties  across  the  axisyrametric 
exhaust  plume.  The  numerical  quadrature  scheme  is  described  in  Appendix 
A.  In  addition,  the  values  of  the  Planck  function  and  the  absorption 
coefficient,  and  their  various  derivatives,  were  evaluated  at  each 
point.  The  equations  used  for  this  purpose  are  given  in  Appendix  B, 


along  with  some  typical  results  to  indicate  the  behavior  of  the 
derivative  terms,  which  define  the  sensitivity  of  the  absorption  to 
temperature  and  species  fluctuations. 

This  type  of  expansion  clarifies  the  contribution  of  each  term  to 
the  time-averaged  radiance  and  suggests  the  order  of  the  importance  of 
the  spatial  correlations  that  are  known  to  be  required  because  of  the 
propagational  nature  of  the  problem.  Identification  of  the  order  of  the 
spatial  correlation  contribution  to  the  radiance  is  a  unique  feature  of 
this  study.  Radiance  from  a  fluctuating  optical  path  with  spatial 
correlations  was  formulated  by  Tan  and  Foster.17  In  fact,  their 
formulation  was  evaluated  exactly  in  closed  form  (without  the  need  for  a 
small  perturbation  expansion),  when  Gaussian  probability  density 
functions  for  the  fluctuations  were  assumed.  However,  they  pointed  out 
that  non-physical  results  (average  transmittance  >  1  at  large  optical 
depths)  were  possible  with  that  assumption.* 


EXAMINATION  OF  THE  CORRELATION  TERMS 

The  precise  contribution  of  the  terms  involving  the  two-point 
correlations  along  the  line  of  sight  requires  the  detailed  numerical 
evaluation  to  be  given  in  the  subsequent  sections.  However,  it  ia 
possible  to  estimate  their  magnitude  and  dependence  on  correlation 
length  and  optical  depth  in  advance.  The  integrals  can  be  evaluated 
analytically  for  the  special  case  where  both  the  mean  and  fluctuating 
properties  are  uniform. 


N  e*p(-  Is  -  s'l/»)as'U-  ij 
m  - 


1  -  exp(-  s/  Z)  dt 


r/Tan,  E. ,  and  P.J.  Foster,  "Radiation  Through  a  Turbulent 
Medium,"  Sixth  International  Heat  Transfer  Conference, 

Toronto,  Canada,  V.  3,  1378,  pp.  403-408.  (paper,  UNCLASSIFIED.) 

•This  nonphysicaily  allowed  result  is  due  to  the  tails  of  the  Gaussian 
distribution  which  admit  small  or  negative  values  of  the  variable  with 
finite  probability,  when  in  reality,  there  is  a  lower  bound  below  which 
the  probability  density  is  zero. 
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=  *<(1  -  e-tL) 


M 


-1  +  kp.A 


[l  -  sxp^- 


“  +  UkQ» 


(19) 


e*“TL)j 

f 

(20) 


The  correlation  coefficient  has  introduced  an  additional  length 
scale  into  the  formication  in  addition  to  the  total  path  length  L,  and 
the  absorption  length,  k^1.  The  nature  of  this  additional  length  is 
perhaps  more  evident  if  we  Took  at  the  limiting  cases  of  3mall  and  large 
optical  paths,  The  asymptotic  oehavior  of  these  integrals  is  given 
in  the  following  table. 


Eq.  (19) 

l  «  L  Jl  »  L 

\  0  1/2  k0L2  (independently  of  Jl) 

tl  +  ®  *  kj1 


Eq.  (20) 

Jl  «  L  Jl  »  L 

1/3  k0L3  (independently  of  Jl) 
2  2k32 


These  limits  suggest  that  for  optically  thin  flows,  (tl  0) ,  the 
contribution  of  the  correlation  terms  is  independent  of  the  correlation 
coefficient  and  depends  only  on  the  total  optical  depth,  i^,  and  total 
geometric  path  length,  L.  In  the  optically  thick  limit  (tl  <*») ,  the 
significance  of  the  additional  length  scale,  Jl,  becomes  apparent.-  When 
J l  is  small,  it  defines  the  order  of  the  contribution  of  the 
correlations.  However,  when  Jl  becomes  large  (compared  to  L),  the 
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absorption  length,  kg1,  defines  the  contribution,  in  agreement  with  the 
well-known  limiting  behavior  of  optically  thick  media.  These  order  of 
magnitude  arguments,  in  particular,  the  absence  of  an  ^dependence  in 
the  optically  thin  limit,  are  verified  in  the  subsequent  numerical 
examples. 

It  should  be  pointed  out  that  we  are  dealing  with  one  form  of  the 
correlation  coefficient  [Eq.  (15)],  and  the  precise  form  of  the 
integrals  and  their  dependence  on  the  correlation  length,  *  are 
particular  to  this  choice.  However,  we  expect  that  the  functional 
dependence  suggested  by  the  asymptotic  limits  in  the  above  table  are 
more  general ,  and  would  hold  for  any  physically  meaningful  correlation 
function . 


RADIATION  FROM  A  TURBULENT  EXHAUST  PLUME 


ISOLATED  SPECTRAL  LINE 

We  anticipate  that  the  predictions  of  the  turbulence  augmented 
radiance  from  the  exhaust  plume  will  depend  on  the  particular  line  of 
sight  through  the  flow  (temperature  and  turbulence  levels),  the  spectral 
location,  the  type  of  line  (transition  between  low  or  high  energy 
level 3) ,  and  the  spectral  resolution.  In  order  to  examine  these  issues 
separately,  we  present  the  radiance  for  an  isolated  spectral  line  before 
dealing  with  the  competing  effects  of  many  lines  in  a  broader  bandpass. 
Typical  exhaust  plume  infrared  signatures  are  usually  low  resolution 
(2-5  cm-1  at  the  highest),  and  do  not  resolve  the  features  of  a  single 
line.  However,  combustion  diagnostics  can  be  high  resolution,  and  the 
results  for  a  single  line  are  of  interest  for  this  application  as  well 
as  illustrating  more  general  features. 

As  a  specific  example,  we  choose  a  broadside  observation  of  the 
plume  at  an  axial  location  40  radii  downstream  from  the  nozzle  exit. 
The  radial  profiles  of  mean  and  fluctuating  properties  at  this  location 
are  shown  in  Figure  14.  This  location  is  in  the  region  of  peak 
afterburning  (Figure  7).  Predictions  of  the  spectral  radiance  from 
isolated  lines  with  different  strengths  and  energy  levels  were  made  for 
lines  of  sight  passing  through  the  middle  of  the  plume  and  also 
displaced  from  the  axis  of  symmetry  (r  =  0.4  rp) .  The  latter  line  of 
sight  emphasizes  the  region  of  increased  turbulence  and  lower 
temperature  that  leads  to  the  largest  effect  of  the  turbulent 
fluctuations.  A  typical  result  at  u>  =  2200  cm-1  (4.5pm)  is  given  in 
Figure  15.  The  line  shapes  for  the  radiance  from  the  mean  properties 
and  from  the  turbulence  augmented  radiance  are  given  in  terms  of  the 
spectral  displacement  Aoo  =  u)  -  u>o  from  the  line  center  normalized  by  the 
half-width,  oq,  at  the  reference  temperature  T0  =  296  K.  Predictions  of 
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x/rj  *40 
«0  *  2220  cm 
tl*5.4 


FIGURE  15.  Spectral  Structure  of  Isolated  Lorentz 
Line.  Broadside  Line  of  Sight  Off  the  Plume  Axis 
at  x/rj  =  40,  A  =  4.03  cm. 
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the  turbulence  augmented  mean  radiance  are  given  for  a  very  large 
correlation  length  (perfect  correlation),  a  correlation  length  equal  to 
the  local  integral  velocity  scale  A,  (our  best  estimate  of  the  true 
correlation  length),  and  a  small  correlation  length  (completely 
uncorrelated  scalar  fluctuations).  These  results  are  for  a  line  of 
moderate  strength  with  a  lower  energy  state  E0  =  4500  cm-1,  typical  of 
lines  that  become  important  at  temperatures  of  600  K  or  more.  In  thi3 
particular  example,  the  turbulence  augmented  radiance  integrated  over 
the  entire  line  width  is  more  than  a  factor  of  four  greater  than  that 
evaluated  at  the  mean  properties.  An  additional  important  result  is 
that  the  integrated  line  radiance  is  not  particularly  sensitive  to  the 
correlation  length,  but  the  precise  line  shape  is  very  strongly 
dependent  on  the  correlation.  This  is  a  result  that  is  potentially 
important  for  any  measurement  that  relies  on  either  the  shape  or  peak 
intensity  of  single  spectral  lines  (i.e.,  line  reversal  temperature 
measurements,  or  rotational  temperature  measurements  from  relative  line 
intensities) . 

A  similar  prediction  for  a  line  of  sight  passing  through  the  plume 
axis  is  given  in  Figure  16.  This  line,  which  is  identical  in  strength 
to  that  used  in  Figure  15,  is  totally  self-absorbed  along  this  line  of 
sight.  Again,  the  integrated  radiance  is  not  significantly  different 
for  the  range  of  correlation  lengths,  but  the  line  structure  is  strongly 
affected.  For  this  choice  of  line  strength  and  lower  energy  level,  the 
absorption  is  intense  and  the  resulting  radiance  is  localized  to  that 
region  of  the  plume  where  the  turbulence  intensity  and  temperature 
combine  to  produce  a  line  structure  that  is  very  different  from  that  of 
a  normal  emission  line,  even  when  it  is  3trongly  self-absorbed. 

Transition  from  completely  correlated  (A  -►  »)  to  totally 
uncorrelated  scalar  fluctuations  is  shown  for  these  two  lines  of  sight 
in  Figure  17.  In  both  these  cases,  the  radiance  from  the  mean 
properties  is  smaller  than  that  with  the  fluctuations  at  all  correlation 
lengths.  This  comparison  also  points  out  that  the  off-axis  line  of 
sight  (r  =  0.4  rp)  through  the  cooler,  more  turbulent  portion  of  the 
flow  has  a  radiance  that  i3  more  intensely  augmented  by  the  turbulent 
fluctuations,  although  its  absolute  level  is  reduced  from  that  of  the 
line  of  sight  through  the  axis.  That  condition  remains  true  for  all 
calculations  performed  in  this  report,  regardless  of  axial  location  or 
spectral  resolution. 

These  predictions  are  for  a  relatively  strong  spectral  line. 
Similar  calculations  for  a  much  weaker  line  show  no  effect  of 
correlation  length,  in  agreement  with  the  limiting  cases  of  the 
correlation  integrals  identified  for  Eqs.  (19)  and  (20).  Moreover,  for 
an  optically  thin  plume,  the  line  of  sight  through  the  axis  shows 
essentially  no  contribution  of  the  turbulent  fluctuations  to  the 
radiance  at  axial  positions  where  the  peak  turbulent  fluctuations  do  not 
coincide  with  the  peak  temperature  (i.e.,  Figure  14).  In  this 
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x/rj  »  40 
w0  =  2220  cm"' 


FIGURE  16.  Spectral  Structure  of  Isolated 
Lorentz  Line.  Broadside  Line  of  Sight 
Through  the  Plume  Axis  at  x/r-s  =  40, 

A  *  4.03  cm. 
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Correlation  Length,  cm 

FIGURE  17.  Integrated  Radiance  of  an  Isolated  Spectral  Line  for 
Various  Correlation  Lengths.  Broadside  Line  of  Sight  at 
x/r.j  =  40.  tp(x)  -  0.429  m. 
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situation,  the  essentially  steady,  but  very  intense  radiance  from  the 
peak  temperature  region  is  affected  very  little  by  the  contributions 
from  the  lower  temperature  outer  region  where  the  fluctuations  are 
relatively  more  intense.  A  line  of  sight  through  the  outer  edge  of  the 
plume  [say  at  r  =  0.4rp(x)]  shows  a  substantial  turbulence  augmentation 
(factor  of  two  for  a  weak  line  with  the  same  lower  energy  level  as  that 
in  Figures  15  and  16)  independently  of  correlation  length. 

An  additional  parameter  that  is  as  important  as  the  optical  depth  i3 
the  lower  energy  level  of  the  particular  spectral  line.  This  affects 
the  temperature  sensitivity  of  the  absorption.  For  example,  the 
Bcltzman  factor  defining  the  population  of  states  with  the  energy  level 
E0  is  exp(-hcE0/KT) .  The  temperature  derivatives  of  the  absorption 
coefficient,  TSk/ST/k,  ^and  T262k/3T2/k  (see  Appendix  B)  therefore 
involve  terms  like  (hcE0/KT)  and  (hcE0/KT)2  and  the  temperature 
sensitivity  is  increased  with  increasing  EQ.  Other  factors,  such  as  the 
line  strength  and  temperature  derivatives  of  the  vibrational  partition 
function  are  essentially  the  same  for  every  line.  Sample  calculations 
of  the  temperature  sensitivity  for  two  typical  energy  levels  and  several 
temperatures  are  given  in  Appendix  B.  These  calculations  suggest  that 
the  temperature  derivatives  change  magnitude  and  sign  according  to  the 
particular  lower  energy  state  and  temperature.  These  effects  are 
somewhat  obscured  by  a  non-homogeneous  path.  However,  the  line  shape 
and  the  effect  of  correlation  length  from  a  line  with  an  energy  level  of 
E0  =  1500  cm-1  are  shown  in  Figures  18  and  19,  respectively,  and  are 
different  from  the  similar  predictions  with  E0  =  4500  cm”1.  This  low 
energy  level  line  (typical  of  those  which  make  up  a  large  portion  of  the 
low  temperature  spectrum  of  CO2  in  the  4.3  pm  region)  is  not 
particularly  sensitive  to  turbulent  fluctuations.  Even  in  the  optically 
thick  region  of  the  line  center  (Figure  18),  the  turbulent  fluctuations 
have  little  effect  on  the  line  shape,  in  contrast  to  the  lines  that 
originate  from  transitions  between  higher  energy  levels  (Figure  15). 

Finally,  we  present  an  indication  of  the  turbulence  effect  on  a  line 
with  a  high  energy  level,  E0  =  7000  cm”1.  The  effect  of  the  correlation 
length  on  the  integrated  line  radiance  is  given  in  Figure  20  for  the  two 
lines  of  sight.  For  this  class  of  lines,  the  temperature  sensitivity  is 
high,  but  the  line  strengths  are  weak,  and  so  the  optical  depths  are 
small.  The  sensitivity  to  the  correlation  length  is  therefore 
diminished.  It  is  interesting  to  note  that  for  this  class  of  lines,  the 
turbulence  augmented  integrated  line  radiance  for  an  uncorrelated  path 
(A  -*■  0)  is  actually  less  than  that  evaluated  at  the  mean  properties. 
However,  the  turbulence  augmented  radiance  for  the  correct  correlation 
length,  A  =  A,  is  greater  than  that  for  the  mean  properties. 

These  few  illustrations  suggest  the  importance  of  the  turbulence 
augmentation  of  the  mean  radiance  from  the  individual  spectral  lines  of 
C02  in  a  small  afterburning  missile  exhaust  plume.  Both  the  line  shape 
and  integrated  radiance  are  affected.  The  extent  of  the  effects  depends 
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FIGURE  18.  Spectral  Structure  of  Isolated 
Lorentz  Line.  Broadside  Line  of  Sight  Through 
the  Plume  Axis  at  x/rj  =  40.  A  =  4.03  cm. 
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FIGURE  19.  Spectral  Structure  of  Isolated 
Lorentz  Line.  Broadside  Line  of  Sight  Off  the 
Plume  Axis  at  x/rj  =  AO,  A  a  A. 03  cm. 
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FIGURE  20.  Integrated  Radiance  of  an  Isolated  Spectral  Line  for 
Various  Correlation  Lengths.  Broadside  Line  of  Sight  at 
x/r.  =  40,  r  (x)  =  .429  m. 


52 


:.4  — 


■I;. 


NWC  TP  6307 


on  optical  depth,  the  particular  line  (energy  level),  and  the 
distribution  of  properties  along  the  optical  path  and  the  extent  of 
correlation  of  these  quantities.  The  largest  contribution  to  the 
augmentation  comes  from  the  additional  fluctuating  source  terms  [term  b 
in  Eq.  (18)].  However,  the  correlations  between  scalar  properties 
(T*,c')  along  the  line  of  sight  can  certainly  be  important.  Past 
studies13*16*17  have  correctly  identified  the  importance  of  the 
fluctuating  source  term;  the  present  3tudy  has  identified  the  explicit 
contribution  of  the  correlations. 


IN-BAND  RADIANCE  -  C02  BLUE  SPIKE 

For  exhaust  plume  signature  applications,  a  spectral  interval  wider 
than  a  single  line  is  of  interest.  The  highest  resolution  band  models 
deal  only  with  0.2  -  2.0  cm-1  resolution18  and  the  most  widely  used 
compilation  is  5  cm-1  resolution. 11  We  chose  a  5  cam-1  bandpass  for  use 
here. 

An  operationally  important  spectral  region  is  the  C02  blue  spike, 
located  between  about  2380  and  2397  cm-1.  This  emission  propagates 
through  long  atmospheric  paths  and  is  a  part  of  the  spectral  signature 
from  all  tactical  missile  exhaust  plume3.  However,  it  is  not  the  most 
sensitive  part  of  the  spectral  signature  to  the  presence  of  turbulent 
fluctuations  because  of  its  reduced  temperature  sensitivity. 16  However, 
because  of  its  technical  importance,  we  have  made  predictions  of  the 
turbulence  augmented  emission  in  the  5  cm-1  bandpass  between  2385  and 
2390  cm"1. 

The  prediction  of  the  emission  from  more  than  a  single  spectral  line 
requires  a  compilation  for  all  the  lines  of  interest  for  the  given 
bandpass.  In  this  region,  we  used  the  AFGL  Line  Atlas.14  Normally,  this 
compilation  is  incomplete  for  high  temperatures.  It  does  not  contain 
the  lines  that  become  important  at  high  temperatures  (E0  5  3000  cm-1). 
However,  these  lines  are  adequate  for  this  particular  spectral  region  up 
to  about  1500  K.  (This  was  verified  by  comparing  the  average  absorption 
computed  with  the  AFGL  Atlas  over  5  cm-1  and  2  cm-1  intervals  with  that 
presented  in  References  11  and  18,  respectively.)  The  use  of  this 
particular  line  atlas  will  underestimate  the  turbulence  augmentation  in 
the  higher  temperature  regions,  but  it  is  adequate  for  the  purpose  of 
this  assessment. 


18Ber stein,  L.S. ,  D.C.  Robertson,  and  J.A.  Conant,  "Band  Model 
Parameters  for  the  4.3  pm  C02  Band  from  200  to  30Q0°K  - 
II.  Predictions,  Comparison  to  Experiment,  and  Application  to  Plume 
Emission  -  Absorption  Calculations,"  J.  Quart.  Spectrosc.  Radiative 
Transfer,  Vol.  23,  1980,  pp.  169-185. 
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Predictions  of  the  tactical  missile  exhaust  plume  radiance  and 
intensity  in  the  2385  -  2390  cm”1  bandpass  were  made  for  various  axial 
locations  in  the  plume.  At  each  axial  location,  the  correlation  length 
was  chosen  to  be  the  local  turbulence  length  scale,  A.  The  length  scale 
and  the  plume  radius  at  the  axial  locations  at  which  these  calculations 
were  made  are  given  below. 


Table  3  -  Axial  Locations 
and  Plume  Dimensions 
for  the  Radiance  Calculation 
r  j  =  .065  m 


x/r  j  A,  m  rp(x) ,  m 


0 

.022 

0.13* 

2 

.024 

.169 

4 

.025 

.182 

10 

.0253 

.247 

20 

.027 

.299 

40 

.0403 

.429 

100 

.093 

.533 

200 

.114 

.585 

400 

.114 

.585 

•Pressure  equilibrated 


Radial  profiles  of  the  mean  and  second-order  correlations  of  the 
fluctuating  properties  have  already  been  given  for  x/rj  =  40  in 
Figure  14.  We  include  additional  profiles  at  the  axial  locations 
x/rj  =  4,  20,  100,  and  200  in  Figures  21,  22,  23.  and  24,  respectively. 
These  profiles  are  given  to  show  the  development  and  spread  of  the  mean 
and  turbulence  quantities  used  in  the  subsequent  radiance  predictions. 

Radial  distributions  of  the  in-band  radiance  in  the  2385  -  2390  cm-1 
bandpass  are  given  at  the  three  axial  locations  x/r^  =  10,  40,  and  100 
in  Figures  25,  26,  and  27,  respectively.  These  locations  are  upstream, 
within,  and  downstream  of  the  peak  afterburning  region.  In  all  cases, 
the  turbulent  fluctuations  have  a  negligible  effect  on  the  mean  radiance 
near  the  plume  axis.  As  mentioned  earlier,  the  radiance  in  the  region 
near  the  axis  is  dominated  by  the  high  temperature,  relatively  less 
turbulent  part  of  the  mixing  region.  The  turbulence  contribution 
becomes  more  important  for  the  off-axis  lines  of  sight  and  the 
downstream  portion  of  the  plume  where  the  peak  turbulence  fluctuations 
and  peak  temperatures  more  nearly  coincide.  For  the  three  axial 
locations  shown,  the  station  radiation  (the  radiation  per  unit  plume 
length)  is  augment? J  by  turbulence  by  the  amounts  of  22%,  20%,  and  37%, 
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FIGURE  23.  Radial  Profiles  of  Mean  and 
Fluctuating  Properties  in  the  Model 
Afterburning  Plume  at  x/rj  =  100.  RSL 
Predictions,  Nominal  Start  Conditions 
and  Turbulence  Length  Scale. 
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FIGURE  24.  Radial  Profiles  of  Mean  and 
Fluctuating  Properties  in  the  Model 
Afterburning  Plume  at  x/rj  =  200.  RSL 
Predictions,  Nominal  Start  Conditions 
and  Turbulence  Length  Scale. 
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FIGURE  27.  Radial  Distribution  of  In-Band  Source 
Radiance  Across  the  Plume  at  x/r^  *  100. 
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at  x/ r j  -  10,  40,  and  100,  respectively. 

The  increased  radiance  ovtr  the  entire  plume  length  is  indicated  by 
the  axial  distribution  of  station  radiation  in  Figure  28.  This  clearly 
shows  the  increasing  relative  importance  of  the  turbulence  on  plume 
emission  with  increasing  axial  distance.  For  this  particular  exhaust 
plume  and  the  2385  -  2390  cm-1  bandpass,  the  turbulence  alimented  mean 
intensity  is  only  24  percent  greater  than  that  calculated  using  only  the 
mean  plume  temperature  and  species  distribution.  This  augmentation  is 
comparable  to  the  uncertainties  usually  associated  with  total  plume 
signature  predictions  under  the  most  favorable  circumstances. 

The  turbulence  augmentation  assumes  a  potentially  more  important 
role  in  a  spatially  resolved  signature.  The  true,  turbulence  augmented 
radiance  is  significantly  larger  than  that  evaluated  at  the  mean 
properties  near  the  outer  edge  of  the  plume.  This  result  was  emphasized 
for  the  single  spectral  lines  discussed  previously,  and  it  is  equally 
true  for  a  wider  bandpass.  The  spectral  radiance  in  the  5  cm-1  bandpass 
is  shown  for  an  off-axis  line  of  sight  (r  =  0.4rp)  at  the  four  axial 
locations  x/rj  =  10..  20,  40,  and  100  in  Figures  29,  30,  31,  and  32, 

respectively.  At  *nese  locations,  the  in-band  radiance  is  increased  by 
72,  45,  100,  and  250  percent,  respectively  over  that  evaluated  without 
the  contributions  of  the  fluctuation.  The  in-band  radiance  is 
substantially  enhanced  for  the  outer  edges  of  this  particular  exhaust 
plume.  Moreover,  we  expect  this  conclusion  to  hold  for  all  afterburning 
plumes.  We  point  out  that  for  plumes  at  higher  altitude  where 
afterburning  is  reduced,  the  entire  plume  is  essentially  a  chemically 
frozen  shear  flow  resembling  the  outer  edge  or  downstream  decaying 
region  of  an  afterburning  plume.  The  turbulence  augmentation  is 
potentially  a  much  larger  contribution  to  the  total  plume  signature  in 
that  case.  We  suggest  that  this  is  the  situation  for  which  Rhodes13 
obtained  such  a  large  effect  of  turbulence  (factor  of  200  increase  in 
station  radiation)  on  the  signature  in  the  2.7  Mm  water  vapor  band. 

Predictions  like  these  were  made  with  both  very  large  and  very  small 
correlation  lengths  in  addition  to  our  best  estimate  of  i  =  A.  The 
results  were  insensitive  to  the  correlation  length.  The  reason  for  this 
is  suggested  by  the  results  for  a  single  spectral  line.  The  correlation 
effect  is  important  in  the  turbulence  augmented  line  structure,  but  not 
in  determining  the  integrated  radiance.  Also,  the  lines  with  lower 
energy  levels  (Figures  18  and  19)  are  less  sensitive  to  the  correlation 
effect.  In  the  2385  -  2390  cm-1  spectral  region,  most  of  the  lines  are 
of  this  class.  The  importance  of  the  correlation  of  fluctuations  along 
the  line  of  sight  is  therefore  not  important  in  this  particular  spectral 
region,  and  the  turbulence  augmentation  is  due  entirely  to  the 
additional  fluctuation  induced  source  terms  [Term  b  in  Eq.  (18)]. 
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FIGURE  28.  Axial  Distribution  of  In-Band  Source 
Station  Radiation  for  the  Model  Afterburning 
Exhaust  Plume.  CO2  Blue  Spike  Region. 
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FIGURE  29.  Source  Spectral  Radiance  for  a  Line  of 
Sight  at  the  Axial  Location  x/rj  =  10  and  Off  the 
Plume  Axis  at  r  =  0.4rp.  C02  Blue  Spike  Region. 
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FIGURE  31.  Source  Spectral  Radiance  for  a  Line 
of  Sight  at  the  Axial  Location  x/rj  =  40  and 
Off  the  Plume  Axis  at  r  =  0.4r  .  CO?  Blue  Spike 
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FIGURE  32.  Source  Spectral  Radiance  for  a  Line  of 
Sight  at  the  Axial  Location  x/rj  =  100  and  Off  the 
Plume  Axis  at  r  *  0.4rp.  CO?  Blue  Spike  Region. 
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Even  for  this  least  sensitive  spectral  region,  the  effect  of  the 
turbulent  fluctuations  is  very  important  on  both  the  magnitude  and 
spectral  detail  of  the  emission  from  the  edges  and  downstream  portion  of 
the  flow.  The  line  shape  and  relative  intensities  of  the  spectral  lines 
are  altered  by  the  turbulence.  These  effects  are  potentially  important 
for  spatially  resolved  signature  applications  and  in  diagnostic 
techniques  relying  on  line  shape  or  relative  intensities. 


IN-BAND  RADIANCE  -  C02  RED  SPIKE 

An  additional  spectral  interval  of  the  C02  4.3  pm  band  of  technical 
interest  is  the  long  wavelength  side  of  the  band  near  4.5  pm 
(2200  cm~l).  This  is  the  spectral  location  of  the  peak  radiance  that  is 
transmitted  through  the  atmosphere.  We  chose  the  5  cm*1  spectral 
interval  between  2220  and  2225  cm”1.  In  this  spectral  region,  the  AFGL 
Atlas  is  inadequate.  Bernstein19  has  shown  that  lines  originating  from 
energy  levels  of  15,000  cm"1  or  more  and  line  densities  of  1000  per  unit 
wavenumber  are  necessary  to  provide  the  correct  absorption  at  high 
temperatures  in  this  spectral  region.  We  utilized  a  synthetic  line 
atlas  for  this  narrow  spectral  region  that  was  developed  at  A.R.A.P.15 
This  formulation  is  more  approximate  than  that  developed  by  Bernstein. 
It  does  not  go  to  as  high  energy  levels  or  contain  as  many  lines. 
However,  the  absorption  averaged  over  2  cm-1  and  5  cm-1  computed  with 
this  synthetic  compilation  compared  favorably  with  Bernstein,  et  el.18 
and  the  NASA  compilation11  up  to  temperatures  of  1500  K.  At  higher 
temperatures,  it  underestimates  the  absorption  and  the  temperature 
sensitivity  of  the  absorption.  Predictions  given  here,  like  those  for 
the  2385  -  2390  cm"1  blue  spike  region,  will  therefore  underestimate  the 
radiance  and  the  sensitivity  of  the  radiance  to  turbulent  fluctuations 
in  the  higher  temperature  regions  of  the  flow. 

Predictions  for  the  2220  -  2225  cm"1  bandpass  were  made  at  the  same 
axi'al  locations  as  those  shown  earlier  for  the  blue  spike  region.  A 
radial  distribution  of  the  in-band  radiance  at  x/rj  =  100  is  given  in 
Figure  33.  Like  the  blue  spike  region,  the  turbulent  fluctuations  have 
no  effect  on  the  mean  radiance  near  the  plume  axis.  This  portion  of  the 
radiance  is  primarily  from  the  high  temperature,  relatively  less 
turbulent  central  portion  of  the  plume.  The  relative  importance  of  the 
turbulent  fluctuations  increases  towards  the  outer  «;dge  of  the  plume, 
but  the  absolute  level  is  sufficiently  reduced  that  the  contribution  of 
the  outer  edges  to  the  station  radiation  is  small.  For  example,  the 
station  radiation  at  the  three  axial  positions  x/rj  =  10,  40,  100  is 

increased  only  by  10,  0,  and  24  percent,  respectively.  These  increases 
are  even  less  than  those  found  for  the  blue  spike  region. 


19.  Bernstein,  L.S.,  "Band  Model  Parameters  for  the  Parallel  Bands  of 
Linear  Triatomic  Molecules  -  I.  Theory,"  J.  Quart.  Spectrosc.  Radiative 
Transfer,  Vol.  23,  1980,  pp.  157-167. 
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FIGURE  33.  Rad>al  Distribution  of  In-Bano 
Source  Radiance  Across  the  Plume  at  x/rj  »  100 
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The  augmentation  of  the  plume  radiance  off  the  axis  of  the  plume  is 
larger  in  the  2220  -  2225  cm-1  bandpass  than  for  the  blue  spike  region. 
For  example,  at  r  =  0.4rp,  the  in-band  radiance  is  increased  by  factors 
of  3.1,  2.4,  and  2.4  for  the  three  axial  locations  x/rj  =  10,  40,  and 
100.  The  augmentation  due  to  turbulence  is  nearly  three  times  greater 
in  the  2220  -  2225  cm**1  bandpass,  as  expected  from  the  greater 
temperature  sensitivity  of  this  region  of  the  spectrum  at  temperature 
levels  that  correspond  to  afterburning  exhaust  plumes.  However,  the 
mean  radiance  also  has  this  sensitivity  and  the  radiance  near  the  plume 
axis  is  sufficiently  large  to  make  the  turbulence  augmentation  a  small 
addition  in  this  particular  wavelength  and  for  the  tactical  missile 
class  of  plumes. 

The  relative  importance  of  the  turbulent  fluctuations  over  the 
entire  plume  length  for  the  2220  -  2225  era"*  bandpass  is  shown  in  Figure 
34.  In  this  particular  example,  the  total  plume  mean  intensity  with  the 
turbulent  fluctuations  is  only  6  percent  greater  than  that  calculated 
without  the  contributions  from  the  fluctuations.  This  rather  small 
increase  is  the  result  of  the  dominant  contribution  by  the  less 
turbulent,  higher  temperature  inner  regions  of  the  plume  (i.e.. 
Figures  14,  22,  and  23). 

At  the  outer  edges,  the  turbulent  fluctuations  provide  an  enhanced 
radiance.  Spectral  distributions  in  the  2220  -  2225  cm"*  bandpass  at 
x/rj  =  10,  20,  40,  and  100  are  given  in  Figures  35,  36,  37,  and  38, 
respectively.  In  these  spectra,  calculated  at  0.01  cm~*  intervals,  both 
the  intensity  and  structure  of  the  lines  are  changed  by  the  turbulence, 
although  the  individual  features  of  the  spectral  lines  are  not  as 
evident  as  in  the  blue  spike  region.  In  this  spectral  region,  there  are 
many  more  lines  contributing  to  the  radiance  (an  order  of  magnitude  or 
more,  depending  on  temperature),  and  the  individual  features  are 
cbscurred.  As  mentioned  previously  the  augmentation  of  the  in-band 
radiance  for  these  off  axis  linss  of  sight  is  greater  for  the 
2220  -  2225  cm"*  spectral  interval.  A3  with  the  2385  -  2390  can"* 
spectral  interval,  the  augmentation  for  lower  temperature, 
non-sfterburning  plumes  would  be  more  like  these  results  for  the  edges 
of  the  plume  and  the  total  signature  would  be  increased  more  than  that 
calculated  here.  However,  we  expset  that  the  slight  increase  in  total 
signature  predicted  here  for  a  particular  tactical  missile  plume  to  be 
generally  true  for  all  low  altitude,  intense  afterburning  plumes.  We 
emphasize  again  that  the  augmentation  at  the  outer  edges,  while  not 
important  for  the  total  signature,  is  potentially  important  in  spatially 
resolved  signatures  and  diagnostic  measurements. 
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FIGURE  34.  Axial  Distribution  of  In-Band  Source 
Station  Radiation  for  the  Model  Afterburning  Exhaust 
Plume.  CO2  Red  Spike  Region. 
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FIGURE  35.  Source  Spectral  Radiance  for  a  Line  of 
Sight  at  the  Axial  Location  x/r^  =  10  and  Off  the 
Plume  Axis  at  r  =  0.4rp.  CO2  Rea  Spike  Region. 


Spectral  Radiance,  w/ 


FIGURE  36.  Source  Spectral  Radiance  for  a  Line 
of  Sight  at  the  Axial  Location  x/rj  *  20  and  Off 
the  Plume  Axis  at  r  »  0.4rp.  CO  9  Red  Spike  Region 


FIGURE  38.  Source  Spectral  Radiance  for  a  Line  of 
Sight  at  the  Axial  Location  x/r^  =  100  and  Off  the 
Plume  Axis  at  r  *  0.4rn.  CO2  Red  Spike  Region. 
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SUMMARY  AND  CONCLUSIONS 


This  report  has  presented  a  set  of  predictions  for  the  mean  radiance 
from  a  turbulent,  afterburning  exhaust  plume  using  a  consistent  set  of 
predictions  for  the  mean  and  turbulence  properties  in  the  flow,  and  a 
precise  formulation  for  the  turbulence  augmented  meai  radiance. 

The  turbulent  flow  properties  were  obtained  from  A.R.A.P.'s  RSL 
code,  which  predicts  the  mean  and  second-order  correlation  of  turbulent 
fluctuations  using  a  second-order  closure  model.  Comparisons  between 
these  predictions  and  those  for  a  two-equation,  turbulence  kinetic 
energy  model  (the  constant  pressure  turbulent  mixing  portion  of  the 
JANNAF  Standard  Plume  Flowfield  Model)  are  in  essential  agreement  for 
constant  density,  low  speed  flows.  Differences  between  these 
predictions  begin  to  occur  for  flow3  with  density  differences  due  either 
to  temperature  or  light  species.  A  major  reason  for  different 
predictions  is  the  different  length  scale  models  used  in  the  two 
calculations,  and  there  is  some  effect  due  to  different  choices  of 
initial  turbulence  levels  and  the  treatment  of  density  fluctuations  in 
variable  density  flows.  For  the  application  at  hand,  a  primary  result 
is  that  the  common  turbulence  quantity  predicted  by  the  two  code3,  the 
turbulence  kinetic  energy,  is  comparable.  An  additional  result  is  that 
the  fluctuations  of  other  scalar  quantities,  such  as  temperature  and 
species  concentrations  develop  differently  from  the  turbulence  kinetic 
energy.  Therefore,  applications  such  as  radiance  predictions  from 
turbulent  flows  require  that  these  additional  scalar  fluctuations  be 
properly  modeled  in  addition  to  the  turbulent  velocity  fluctuations. 

For  the  particular  flow  examined  here,  and  with  the  global  model  for 
the  afterburning  reaction,  we  found  that  the  effect  of  the  turbulent 
fluctuations  on  the  chemical  reactions  has  only  a  small  influence  on  the 
temperature  and  species  concentrations  in  the  afterburning  region. 

A  formulation  for  the  mean  spectral  radiance  accounting  for  the 
additional  fluctuating  source  and  transmittance  terms  created  by 
fluctuations  in  temperature  and  chemical  species  was  derived.  It  is 
essentially  a  direct  expansion  of  the  monochromic  radiance,  up  to  and 
including  second-order  correlations  between  the  fluctuations.  It  is 
uniformly  valid  for  all  optical  depths  but  requires  that  the  turbulent 
fluctuation  remain  small  compared  to  their  mean  values,  a  condition  that 
is  satisfied  by  the  flowfield  used  here.  A  unique  feature  is  that  the 
contributions  of  spatial  correlations  between  fluctuating  properties 
along  the  line  of  sight  are  explicitly  obtained. 

The  spectrally  resolved  mean  radiance  with  and  without  the 
contributions  from  the  fluctuations  was  computed  for  two  regions  in  the 
4.3  pm  COg  vibrational-rotational  band  for  lines  with  a  Lorentz  line 
shape.  The  spectral  information  for  the  line  strengths  and  positions 
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was  taken  from  the  AFGL  compilation  in  the  2385  -  2390  cm**1  spectral 
region  and  from  a  synthetic  compilation  in  the  2220  -  2225  cm-1 
bandpass.  Both  compilations  were  verified  to  be  adequate  up  to  about 
1500  K,  but  underpredict  the  absorption  at  higher  temepratures. 

Predictions  for  isolated  spectral  lines  were  used  to  determine  the 
importance  of  the  correlation  between  fluctuating  properties  along  the 
line  of  sight.  For  optically  thin  paths,  there  is  no  effect  of 
correlation.  However,  at  intermediate  and  large  optical  depths,  the 
correlations  affect  the  line  structure,  but  do  not  strongly  affect  the 
integrated  line  radiance.  An  important  parameter  governing  the 
sensitivity  of  the  line  to  the  turbulent  fluctuations  is  the  lower 
energy  state  of  the  particular  line.  The  sensitivity  increases  with  an 
increasing  lower  energy  level,  suggesting  that  the  least  sensitive 
absorption  is  due  to  the  lines  that  determine  the  low  temperature 
spectra,  and  that  the  most  sensitive  lines  are  those  which  contribute  at 
the  higher  temperatures.  These  observations  have  implications  to  plume 
diagnostics  that  utilize  properties  of  the  individual  line  shape  or 
relative  intensities  of  single  lines. 

In-band  radiance  for  wider  bandpasses  (5  cm-1),  wa3  predicted  for 
the  C02  blue  spike  (2385  -  2390  cm-1)  and  red  spike  (2220  -  2225  cm-1) 
spectral  regions.  For  the  blue  spike  region,  the  station  radiation  in 
the  afterburning  and  downstream  region  was  increased  by  at  most  37 
percent,  and  the  total  plume  intensity  for  broadside  observation  was 
increased  by  24  percent  by  the  turbulent  fluctuations.  The  small 
augmentation  is  because  the  predominant  contribution  to  the  signature 
comes  from  those  lines  of  sight  passing  near  the  plume  axis,  for  which 
the  turbulent  augmentation  is  very  small.  The  outer  edge3  of  the  plume 
are  substantially  augmented  by  the  turbulence,  by  factors  as  much  as  2.5 
times  the  radiance  without  the  fluctuations.  These  outer  regions  do  not 
contribute  significantly  to  the  total  plume  intensity,  but  could  be 
important  in  a  spatially  resolved  signature. 

In-band  radiance  in  the  red  spike  region,  2220  -  2225  cm-1,  has 
essentially  the  same  behavior  as  that  in  the  blue  spike,  except  that  the 
center  portion  of  the  plume  is  even  more  dominant  and  the  augmentation 
by  the  fluctuations  to  the  total  plume  intensity  is  even  less;  only  6 
percent  for  the  example  used  here.  For  this  bandpass  the  station 
radiation  is  increased  by  at  most  33  percent,  but  only  downstream,  where 
the  effect  on  the  total  intensity  is  small.  The  augmentation  of  the 
radiance  at  the  edges  of  the  plume  is  substantial,  but  not  important  for 
the  total  plume  intensity. 

The  effect  of  the  spatial  correlations  of  the  temperature  and 
species  concentrations  was  very  small  for  the  radiance  in  these  two 
bandpasses.  The  correlations  affect  the  line  structure  but  not  the  line 
strength  and,  therefore,  a  broader  bandwidth  radiance  is  much  less 
sensitive  to  the  correlation.  The  turbulence  augmentation  of  the  wider 
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bandpass  radiance  is  due  primarily  to  the  additionl  radiance  source 
terms  created  by  the  single  point  correlations  of  turbulent 
fluctuations,  and  not  to  the  fluctuation  augmented  transmittance,  to 
which  the  two-point  spatial  correlations  contribute. 

We  conclude  that,  for  this  specific  case  of  a  typical  low  altitude, 
afterburning  exhaust  plume  from  a  tactical  missile,  turbulent 
fluctuations  are  not  important  contributors  to  the  total  plume  intensity 
in  either  the  COg  blue  spike  or  red  spike  spectral  regions.  The 
predicted  augmentations  in  these  spectral  regions  of  24  and  6  percent, 
respectively,  are  well  within  the  uncertainty  of  contemporary  flowfield 
and  signature  predictions.  We  expect  this  conclusion  to  hold  for  the 
entire  class  of  tactical  missile  plumes  at  low  altitude. 

However,  the  turbulence  augmentation  could  be  substantial  for 
non-afterburning  plumes  such  as  higher  altitude  rocket  plumes  and  for 
aircraft  plumes.  For  this  class  of  plumes,  the  flow  is  more  like  the 
outer  edges  and  downstream  portions  of  an  afterburning  plume,  for  which 
the  predicted  turbulence  augmentation  was  a  factor  of  from  two  to  three. 
Turbulence  could  then  substantially  increase  the  total  plume  intensity 
from  these  types  of  exhaust  plumes.  For  the  same  reasons,  turbulent 
augmentation  of  radiation  observed  in  the  forward  hemisphere  is 
expected.  In  cases  of  forward  hemisphere  observation,  much  of  the 
afterburning  region  of  the  plume  is  blocked  by  the  missile  structure, 
and  only  the  outer  edges  and  downstream  portions  contribute  to  the 
observed  signature. 


RECOMMENDATIONS 


The  results  of  this  study  suggest  several  extensions.  There  are 
clear  implications  of  these  results  for  high  spectral  resolution 
diagnostic  techniques  such  as  temperature  measurements  by  line  reversal 
techniques  and  rotational  temperature  measurements  from  the  distribution 
of  line  intensities.  It  would  be  worthwhile  to  examine  the  results  of 
these  types  of  measurements  in  experiments  where  the  flow  is  turbulent 
with  the  turbulent  augmentation  included. 

An  important  contemporary  plume  diagnostic  that  would  be  affected  by 
the  predicted  turbulence  augmentation  in  the  wider  bandpass  (~  5  cm"1) 
is  the  emission/ absorption  inversion  technique  used  to  infer  temperature 
and  radiating  species  profiles  from  axisymmetrie  flows.  In  these 
measurements,  the  entire  transverse  profiles  of  emission  and 
transmission  are  utilized.  The  substantial  turbulence  augmentations  at 
the  outer  edges  of  the  plume  suggest  that  the  actual  inferred  properties 
are  radiation  weighted  mean  temperature  and  species  concentrations  that 
could  differ  from  the  presumed  normal  mean  properties.  An  examination 
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of  experimental  measurements  with  some  indication  of  the  augmentation  of 
the  emission  might  indicate  the  effect  on  the  inferred  temperature. 
This  application  would  require  a  prediction  of  the  time-averaged 
fluctuation  augmented  transmittance,  which  was  not  given  here 
explicitly,  but  could  be  easily  derived. 

Finally,  we  argue  that  it  is  probably  feasible  to  include  a 
prediction  of  the  turbulence  augmented  radiance  in  the  normal  signature 
predictions.  The  missing  element  is  a  formulation  of  the  turbulence 
augmentation  for  band  models.  Rhodes's  ensemble  averaging  approach  is 
probably  too  time  consuming  for  most  applications. 13  This  capability 
would  be  most  valuable  for  higher  altitude  rocket  plumes  or  aircraft 
plumes.  For  line-by-line  predictions,  the  formulation  presented  here  is 
directly  applicable.  In  all  these  applications,  however,  the  primary 
requirement  is  a  plume  flowfield  prediction  that  includes  reliable 
models  for  the  scalar  fluctuating  quantities  that  affect  the  plume 
emission . 
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Appendix  A 

NUMERICAL  EVALUATION  OF  INTEGRALS 

There  are  t vo  types  of  integrals  to  be  evaluated  in  the  expression 
for  the  time-averaged  radianoe.  One  is  the  integral  over  the  optical 
path,  and  the  second  is  the  integral  of  tile  correlation  coefficient  over 
the  geometric  distance  along  the  line  of  3ight.  The  two  integrals  are 
treated  similarly,  but  the  nianerioal  quadratures  are  different.  Each 
will  be  described  separately. 

INTEGRAL  OVER  THE  PATH 

This  integral  is  of  the  form 


(A-l) 


where  F(s)  is  any  form  of  the  integrand  [even  including  other  integrals, 
i.e.,  Eq.  (18)].  The  independent  variable  is  transformed  to  optical 
depth 


I 


■/ 


x  =  I  k  cds' 
o 


J> 


F(x)  exp(-  x)dx 


(A- 2) 


(A- 3) 


The  optical  path  can  vary  between  very  small  and  very  large  values, 
0  <  x  <  ®.  Because  of  the  exponential  weighting,  a  single  quadrature 
scheme  may  not  be  uniformly  accurate.  For  example,  if  xl  «  1,  then  the 
entire  optical  path  makes  a  contribution  to  the  integral,  while  if 
xj,  »  1 ,  only  che  thin  layer  of  plume  closest  to  the  observer  determines 
the  radiance.  In  order  to  correctly  account  for  these  effects  in  a 
uniformly  accurate  manner,  the  numerical  evaluation  of  Eq.  (A-3)  was 
performed  with  two  quadrature  schemes,  depending  on  the  value  of  the 
total  optical  depth. 
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Limits  of  Integration  Integration  Scheme 


\  <  10 

\  >  10 


0  <  t  < 


0  <  x  <  « 


Lobatto 


Laguerre 


In  the  case  where  x^  <  10,  the  variable  of  integration  is  5  5  t/t l»  an<* 
the  integral  is 


■■jP 


F(x)  e*“t  dx  =  £  F(t^)exp(-5jx^)wi 


(A-4) 


where  ^  are  uke  roots,  and  w^  the  weights  of  the  integration, 
large  optical  path  limit, 


For  the 


F(x)  e--1  dx  ~  53  F(0  w^  • 

i=1 


(A-5) 


In  both  cases,  N  =  16  points  [half  of  a  30  point  scheme  for  the  finite 
path  length  case,  Eq,  (A4)J.  These  numerical  quadratures  were  highly 
accurate  ahd  gave  essentially  enact  results  for  all  optical  depths  when 
compared  to  cases  that  could  be  evaluated  analytically  [constant  values 
for  the  integrand  F(s)].  The  results  were  insensitive  to  the  choice  of 
cut-off  between  the  finite  ar.d  infir  Ite  optical  depth  limits. 

INTEGRALS  OVER  THE  SPATIAL  CORRELATION  COEFFICIENT 

These  integrals  are  of  the  form 


F(s*>  exp(-  Is  -  s’  I / A) d s * 


(A~6) 


Again,  because  of  the  exponential  weighting,  an  arbitrary  numerical 
quadrature  (such  as  a  trapezoidal  rule)  is  not  uniformly  accurate  unless 
the  increments  are  always  3mall  compared  to  i,  the  correlation  length. 
To  circumvent  this  inconvenience,  the  variable  of  integration  was 
tr an formed  to 

l  =  exp(-  Is  -  s' i /A)  ,so  that 
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I  * 


F(£)  d£ 


(A-7) 


This  integral  was  evaluated  by  a  normal,  trapezoidal  rule.  The 
numerical  evaluation  gave  essentially  exact  results  for  arbitrary  values 
of  SL.  It  is  this  convenient  transformation  that  mitigated  against  us 
examining  other  forms  of  correlation  coefficients. 
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Appendix  B 

DERIVATIVES  OF  THE  PLANCK  FUNCTION  AND 
SPECTRAL  ABSORPTION  COEFFICIENT 


PLANCK  FUNCTION 


C^3 

B  (T)  =  - - - 

exp(hco/KT)  -  1 


(B-l) 


Let  X  =  hco/KT,  then 


T  dB/dT  =  BXeX/(exp  X  -  1) 


(B-2) 


and 


T2  d2B/dT2  s  2(T  dB/dT) (T/B  dB/dT  -  1  -  X/2)  (B-3) 

ABSORPTION  COEFFICIENT 

The  spectral  absorption  coefficient  for  a  rotational  line  is 


k  =  S  (T)f(a  ) 


(B-4) 


where  SjCT)  is  the  line  strength,  f(a)  is  the  line  shape  function,  and 
at  the  half-width  for  the  ith  spectral  line  with  center  wavenumber  at 
u>oi.  We  assume  a  Lo^entz  line  shape 


f(tti)  =  1/n 


_ ai _ 

ai  +  (u)  -  uoi) 2 


(B-5) 


The  line  strength  is  given  by 
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S  exp(~  hcfc  /KT ) X 1  -  expt-  hC(j/KT)]Q  (T  )T  /T 
Sjd)  =  pN  -2i - - 2 - v__o_o_  (b_6) 

Qr(T)Qv(T)  exp(-  hSE^/KT0)[1  -  exp(-  h6'(1/KT0)3 


Soi  is  the  line  strength  at  the  reference  temperature  T0,  and  E0  is  the 
energy  of  the  lower  level  of  the  particular  vibrational-rotational 
transition  which  is  responsible  for  the  line.  and  Qv  are  the 
rotational  and  vibrational  partition  functions,  respectively.  The  first 
exponential  term  accounts  for  the  distribution  of  molecules  in  the 
energy  state  E0  (the  Boltzman  factor),  and  the  second  exponential  term 
accounts  for  stimulated  emission.  The  temperature  ratio  T0/T  i3 
explicitly  included  here  to  account  for  the  fact  that  SQ  is  defined  on  a 
per  molecule  basis.*  The  temperature  ratio  is  conventionally  carried  as 
a  separate  term. 

The  partial  derivatives  of  the  absorption  coefficient  are  calculated 
in  the  following  manner. 


ak/aT  =  as/aTf+  saf/aai  aai/aT 


ak/  Sc  : 

■  Saf/5a^  5oti/5c 

a2k 

af 

5a.- 
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a2f /do^N2 

„  df 
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5T2 

dT 
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5T 

5a|\  Si/ 
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d2k 

as  af 
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(B-8) 


(B-9) 

(B-10) 

(B-ll) 


Each  of  those  terras  is  defined  below: 


«k(cm-i)  -  Jc(molecule“l)  •  N  s  kN0(N/No)  =  kN0pT0/T  where  N  is  the 
molecular  number  density  and  p  the  pressure.  We  define  the  absorption 
ooeffeient  to  be  kW0pT0/T. 
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f(otj)  *  otj/i^a?  ♦ 

(B-12) 

5f/dc|  s  f ( 1  —  2i»^f)/ot£ 

(B-13) 

52f/'5a^  =  2«f2(4Ka^f  -  3)/a£ 

(B-14) 

«i  =  «0p(T0/T) 1/2(yiO  *  Yj°j^ 

(B-15) 

T  aa^/ftT  =  -otj/2 

(B-16) 

T2  &2ai/aT2  =  3«i/4 

(B-17) 

dai/ac  =  (jQpdQ/T) 1/2  (y^  -  yj) 

(B-18) 

a^/aodT  =  -i/2  dai/ac 

(B-19) 

Yi  and  y}  are  constant  coefficients  which  account  for  self  and 
gas  broadening.  In  this  work,  yj  =  1.4,  yj  =  1.07,  and  cj  =  1  - 
these  conditions. 

foreign 
c.  For 

d2otj/ac2  =  o 

(B-20) 

Derivatives  of  the  temperature  dependent  portion  of  the  absorption 
coefficient  are  more  complicated.  First,  the  partition  functions  are 
(for  C02): 


rotation;  Qr(T)  =  T/T0 


(B-21) 


vibration: 


QV(T) 


_ 1 _ 

(1  -  exp  -  hJc?U'j/KT)2(l  -  exp  -  ti&ug/KT)(l  -  exp  -  hcf^/KT) 


(B-22) 
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where  wj  =  667.379  cm"1,  w?  =  1388.187  cm**1,  W3  a  2349.146  cm"1. 

Normally,  the  stimulated  emission  term  and  Qv  at  the  reference 
temperature,  T0,  (296°K),  which  are  nearly  unity,  are  not  explicitly 
included  in  Eq.  (86). 

For  the  purpose  of  evaluating  the  derivative  of  S(T),  we  define  the 
following  terms: 


XQ  =  h<5E0/KT 
X  =  hSVKT 
Xt  =  hCty/KT 
X2  =  hSug/KT 
X3  =  her^/KT 


y0  =  X0e~X0/(i 
y  =  Xe“x/(1  - 
yi  =  X  -j  e"*x  1/(1 
y2  =  X2e-*x2/(1 
y3  =  X3e"x3/(1 


-  exo> 
e“x) 

-  e"x1 ) 

-  e“x2) 

-  e"x3> 


then , 


T&S 

dT 


=  S(T)(Xq  -  2  -  y  -  2yt  -  y2  -  y^ 


(B-23) 


T2  ~  s  -T—  +  T— (X  -  2  -  y  -  2y  -  y  -  y  )  • 
aT2  QT  9T  0  12  3 


.  -Xq  +  y(1  -  X  -  y)  +  2y1(1  -  X,  -  y,)  +  yg(1  -  X2  -  y2>  + 


+  y3(1  '  X3  ~  V 


(B-24) 


These  derivatives  determine  the  sensitivity  of  the  absorption  to 
temperature  and  radiating  species  concentration  fluctuations.  The 
primary  parameters  ar?  temperature  and  wavelength,  but  another  important 
parameter  is  the  energy  of  the  lower  state  of  the  transition,  E0.  The 
ratio,  hSE0/KT,  affects  the  temperature  dependence  of  the  absorption. 
It  is  possible  for  lines  with  different  lower  state  energies  to  have 
opposite  temperature  sensitivities.  Some  examples  are  given  in  Figures 
B1  -  B6  for  typical  temperatures  found  in  the  plume  and  for  wave-lengths 
near  4.5  pm.  The  plots  are  normalized  by  the  absorption  at  the  line 
center,  kQ,  and  the  wavenumber  measured  from  the  line  center. 
Aw  =  w  -  uq  is  normalized  by  the  half-width,  ai  =  a.  evaluated  at  the 
temperature,  T. 
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FIGURE  B-l.  Normalized  Derivatives  of  the 
Spectral  Absorption  Coefficient:  Single 
Line,  Lorentz  Line  Shape,  E0  «*  1500  cm-1, 

T  =  300  K. 
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FIGURE  B-2.  Normalized  Derivatives  of  the 
Spectral  Absorption  Coefficient:  Single  Line, 
Lorentz  Line  Shape,  Eq  *  1500  cm-1,  T  =  1000  K. 


r- 
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FIGURE  B-5.  Normalized  Derivatives  of  the 
Spectral  Absorption  Coefficient:  Single 
Line,  Lorentz  Line  Shape,  Eg  =  4500  cm-1, 

T  =  1000  K. 


FIGURE  B-6.  Normalized  Derivatives  of  the 
Spectral  Absorption  Coefficient:  Single 
Line,  Lorentz  Line  Shape,  E0  =  4500  cm"1, 

T  =  2000  K, 
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NOMENCLATURE 

Planck  function,  B0  =  BU(T) 

Scale  constant  for  turbulent  length  scale 

Constant  in  Planck  function, 

C1  =  1.19088  x  10-12  u>-cm2/str 

Species  mole  fraction  (primarily  for  the  infrared 
active  species  CO2) 

Speed  of  light 

Activation  energy 

E/RT 

Lower  energy  level  of  transition  for  an  individual 
spectral  line 

Arbitrary  function 

Lorentz  line  shape  function 

enthalpy 

Planck’s  constant 

Boltzman's  constant 

Forward  reaction  rate  coefficient 

Spectral  absorption  coefficient,  k0  =  k^CT), 
also  denoted  by  k 

Geometric  length  of  optical  path  through  plume 
Correlation  length  for  scalar  fluctuations 
Number  density  of  molecules 
Spectral  radiance 
Pressure 


95 


NWC  TP  6307 


^v.Qr 

q2 

R 

r 

rj 

Vx) 

rTT.  rcT»  rcc 

S 

s 

T 

u 

u' ,  V* ,  w' 
wi 

X 


Vibrational  and  rotational  partition  functions 

Mean  square  turbulence  velocity 

Universal  gas  constant 

Radial  coordinate 

Exhaust  nozzle  radius 

Plume  radius  at  axial  location,  x 

correlation  coefficients  between  fluctuating 
temperature  and  species 

Line  strength 

Distance  along  line  of  sight 

Temperature 

Mean  axial  velocity 

Turbulence  velocities  in  x,  r,  9  direction 
Weights  for  Gaussian  Quadrature 
Axial  coordinate 


a  Reactant  mass  fraction 

oci  Spectral  line  half-width 

p  Reactant  mass  fraction 

Y  Product  mass  fraction 

Yi  Collision  broadening  coefficient 

e  Turbulence  energy  dissipation 

A  Turbulence  velocity  integral  scale 

x  Optical  depth 

p  Fluid  density 
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Wavenumber 

Variable  of  integration 


Subscripts 
jet,  j 

00 

0 

pk 

m 

I  s’  « I  s” 
i.j 


Jet  exit 
Ambient 

Reference  condition  or  evaluated  at  mean  temperature 
Peak  value 
Axis  value 

Evaluated  at  the  location,  s’,  s" 

Chemical  species  i.j 


Superscripts 


Mean  (time-average) 

Fluctuating  (note  the  prime  is  also  used  to  denote 
a  dummy  variable  of  integration 
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